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ABSTRACT 

The electroweak phase transition in the magnetic and hypermagnetic fields is 
studied in the Standard Model on the base of investigation of symmetry behaviour 
within the consistent effective potential of the scalar and magnetic fields at finite 
temperature. It includes the one-loop and daisy diagram contributions. All dis- 
covered fundamental fermions and bosons are taken into consideration with their 
actual masses. The Higgs boson mass is chosen to be in the energy interval 75 
GeV < rriH < 115 GeV. The effective potential calculated is real at sufficiently 
high temperatures due to mutual cancellation of the imaginary terms entering 
the one-loop and the daisy diagram parts. Symmetry behaviour shows that nei- 
ther the magnetic nor the hypermagnetic field does not produce the sufficiently 
strong first order phase transition. For the field strengths H, Hy > 10 23 G the 
electroweak phase transition is of second order at all. Therefore, baryogenesis 
does not survive in the Standard Model in smooth magnetic fields. The problems 
on generation of the fields at high temperature and their stabilization are also 
discussed in a consistent way. In particular, it is determined that the nonabelian 
component of the magnetic field (gH) 1 ^ 2 ~ g 4 ^ 3 T has to be produced sponta- 
neously. To investigate the stability problem the W^-boson mass operator in the 
magnetic field at high temperature is calculated in one-loop approximation. The 
comparison with results obtained in other approaches is done. 



1. INTRODUCTION 

Among interesting problems of nowadays high energy physics there are two ones 
which, at first glance, are not connected with each other. These are the value of 
the Higgs boson mass tuh and the strengths of magnetic fields H which could 
be present in the early universe (see surveys |jj - |§). Both problems are of 
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paramount importance for particle physics and cosmology. For instance, as was 
shown in Refs. |J, ||, |J a large scale homogeneous hypercharge magnetic field 
Hy must essentially influence the type of the electroweak (EW) phase transition 
making it strong first order. Interest to the effects of strong magnetic fields has 
considerably increased recently when it was realized that a standard baryogen- 
esis scenario in the Standard Model (SM) of elementary particles could not be 
established without the fields fl3| . Different mechanisms to produce the fields 
have been proposed @, Jl7|, 0, jn|, |§, &H, 0- One of the purposes 



of the present paper is to discuss some recent results on these topics. At the 
same time, in the Monte Carlo simulations and by the nonperturbative methods 
of quantum field theory |L5] , [14[] the magnetic mass of nonabelian gauge fields of 



order m mag , ~ g 2 T has been determined. It screens the nonabelian component of 
the magnetic field at long distances. This is one of the reasons why an interest to 
the hypermagnetic field Hy was excited. This field, owing to its abelian nature, 
is not screened at finite temperature. The influence of hypermagnetic fields on 
symmetry behaviour has been first investigated for many years ago in Ref. JTE 



where the similarity to superconductivity was emphasized. This important ob- 
servation plays a decisive role in the description of the EW phase transition. The 
problem which requirs further investigations is the generation of strong abelian 
fields at high temperature. Interesting mechanism connecting with an abilian 
anomaly was suggested in Ref. [0. 

Much more involving is situation with the magnetic field, although it is studied 
for many years. Currently, there are not common opinions about either the 
ways of producing and stabilization the field or its role at high temperature (see 
papers [Tj|, P0|, |Z71 and references therein). This situation, probably, finds an 
explanation in the nonabelian nature of the field and lack of simple analogies in 
condensed matter physics (mainly superconductivity is discussed) . As an example 
of peculiarities, let us mention the phenomenon of condensation of the W- and 
Z-boson fields which originates from the instability of the perturbative vacuum in 
strong magnetic fields (see Refs. |JD| - P3). In such a situation, to treat problems 
with the external fields correctly the consistent calculations have to be carried 
out. Recent investigations of the phase transitions in the magnetic fields at high 
temperature ||34jl , p5| , |27|j have used as a qualitative picture of the phenomenon 



the description in Refs. fl36| , |20|] derived at a classical level for the case run = Tnz 
corresponding to second type superconductivity. However, in these papers the 
effects of fermions (light and heavy) as well as the radiation corrections in the 
fields at high temperature have not been included but play an important role. 
The second goal of the present paper is to elaborate the situation with magnetic 
fields at high temperatures. We shall consider all the problems mentioned within 
consistent calculations allowing for the one-loop effective potential (EP) and the 
daisy diagram contributions in the external fields at high temperature. To find 
the latter ones the one-loop polarization functions of gauge bosons under these 
external conditions will be computed. As it is occured, the influence of the fields 
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crucially depends on the values of the particle masses. So, to have a reliable 
quantitative description of the EW phase transition we fix them to be equal to 
the present day experimental data. The Higgs boson mass will be taken in the 
energy interval 75 GeV < mn < 115 GeV. The low bound corresponds to the 
mass values when perturbative methods are reliable. The upper bound is chosen 
to fit the present experimental low limit m# > 90 GeV. In our calculations, the 
external fields will taken into account exactly through the Green functions. This 
insures, the results to be obtained should reproduce correctly the effects of the 
fields for all the values niu considered. It is also important to notice that in 
strong fields at high temperatures light fermions dominate, as it follows from 
the term H 2 logT /rrif entering the one- loop effective potential, where my is the 
fermion mass. Actually, at different temperatures the different fermions give the 
dominant contributions. 

The concept of symmetry restoration at high temperature has been intensively 
used in studying of the evolution of the universe at its early stages. Nowadays it 
is a corner stone in investigations of various problems of cosmology and particle 



physics |37|], |38|. In particular, the type of the EW phase transition and, hence, 
a further evolution of the universe depends on the mass of Higgs boson. As 
we mentioned above, the idea that strong fields make the EW phase transition 
strong first order, that is necessary to retain the standard baryogenesis (see recent 



survey |13]|), requirs further investigations with radiation corrections allowed for. 
Since all the masses of fundamental particles, except m#, are known one is able 
to investigate in detail the phase transition as the function of this parameter and 
to determine the properties of the vacuum. That is the main goal of the present 
paper. 

Various aspects of the phase transitions in magnetic fields at high temperature 
have been discussed in literature - J42|, [[18]. In Refs. f43|, p4| , considering 



the boson part of the Salam- Weinberg model, the EW phase transition in the 
strong fields was investigated and the vacuum structures of the phases have also 
been described. In Ref. [[TJj side by side with temperature and magnetic field 
a chemical potential was incorporated. But the role of fermions has not been 
studied in detail. 

Another aspect of the EW phase transition, which also has not been elabo- 
rated, is the influence of so-called daisy (or ring) diagrams at high temperature 
and strong fields. At zero field it has been investigated in Refs. |4^], |j5], |{46 



where the importance of these diagrams to correctly describe symmetry behaviour 
was emphasized. In Ref. [^6[ the t-quark mass was chosen of order 110 GeV. 
So, to account of the experimental value m t = 175 GeV, it has to be revised. 
In the present paper the EW phase transition will be studied for the case of the 
constant fields Hy and H . This is an adequate approximation for strong fields 
in the cases of the second order phase transition and the initial stage of the first 
order one when the bubbles are not large ||. 

The content is as follows. In sect. 2, for convenience of readers, the results of 
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the investigations announced are described in a qualitative manner. In sects. 3, 
4 the contributions of bosons and fermions to the one-loop EP (T, H, C ) of 
classical the scalar <p c and the external magnetic fields are calculated in the form 
convenient for numeric investigations. In sect. 5 the correlation corrections due 
to daisy diagrams are computed and the vacuum stability condition at high tem- 
perature is discussed. Special attention is devoted to computation of the daisy 
diagrams with the unstable (tachyonic) mode presenting in the magnetic field 
in the VF-boson spectrum. In sect. 6 the restored phase with external fields is 
described. In sect. 7 the high temperature expansion of the EP is present. Then, 
in sect. 8 the EW phase transition in the hypercharge magnetic field is investi- 
gated. The same for the magnetic field case is carried out in sect. 9. In these 
two sections the detailed analysis of the phase transition is done. The one-loop 
polarization functions of M^-bosons at high temperatures and strong magnetic 
fields are calculated in sect. 10. This, in particular, gives possibility to study 
self-consistently the vacuum stability. The spontaneous vacuum magnetization at 
high temperature is investigated in sect. 11. It is found that the magnetic fields 
of order (gH) 1 / 2 ~ g 4 ^ 3 T is generated. Such strong fields affect all the processes 
at high temperatures. Comparison of the results obtained with that of other 
approaches is done in sect. 12. Discussion and concluding remarks are given in 
sect. 13. Appendix contains necessary information on the Mellin transformation 
technique used in calculations of the high temperature asymptotics of the EP. 



2 . QUALITATIVE PICTURE OF THE EW PHASE TRAN- 
SITION IN MAGNETIC FIELDS 

As it is belived nowadays, the presence of different kind strong magnetic fields 
in the early universe is rather resonable than exotic phenomenon. In literature 
on this topic a lot of dynamic mechanisms to generate the fields are proposed 
(see surveys ||, H and references therein). In the present paper, we are not 
going to consider all of them. Our goal is to describe a consistent picture of the 
influence of the fields on the EW phase transition. In this section, we consider in 
a qualitative manner the mechanisms of producing the external fields, the ways of 
their action on the vacuum at zero and finite temperature and present the main 
results of our investigation. 



2.1 The generation of the primordial magnetic fields 



The generation of the fields can be devided in two classes: 1) generation due to 
processes that had happened at the EW phase transition; 2) creation of the fields 
before the EW phase transition epoch (for instance, at a GUT scale). To the 
first class we refer: the rotation of bubbles at the first order phase transition |f22 
fluctuation of gradients of the Higgs field [||; bubble collisions |23 



251 . These 
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mechanisms produce microscopic fields which had to be amplified by magnetohy- 
drodynamic processes up to the macroscopic values required to fit the astronomic 



observations, as it is discussed in Refs. [50], [51 



To the second class of the processes we refer the spontaneous vacuum magne- 
tization at high temperature discovered first at zero temperature by Savvidy |52 



and the fields generated by strings [48| . Ones had been created, these fields were 
frozen in a cosmic plasma, evolved with it during the expansion of the universe 
and present at the EW phase transition. In what follows, the former mechanism 
will be investigated in detail. We shall show that it does work at high tem- 
perature. So, it could serve to produce the seed magnetic field (in contrast to 
conclusions of Refs. @], fl35f). 



Let us remind the results of Refs. p2[, [O], |20[, [O]. The spontaneous vac- 



uum magnetization has been derived from the investigation of the EP of covari- 
antly constant (sourceless) chromomagnetic field H a = H5 a3 which is a solution 
to the classical Yang-Mills field equations, where H = const and a is an isotopic 
index, 

V(H,T) = ^- + V^(H,T). (1) 

It includes the tree-level and the one-loop parts. The minimum of the EP at high 
temperature T corresponds to the nonzero magnetic field of order (gH)^ 1 ^ 2 ^ ~ 
g 2 T, g is gauge coupling constant. [] In the EW theory, the a = 3 component of 
the weak isospin just corresponds to the nonabelian part of usual magnetic field 
which we observe in the broken phase. 

Very important for our analysis is the value of the vacuum field strength 
~ g 2 T. In fact, as it will be shown in sect. 11, this value is increased when 
the correlation corrections are taken into acccunt. The field is screened at long 
distances / > l/m mag ~ l/g 2 T by the magnetic mass of gauge field. However, 
inside this space domain the strong fields may exist and affect all the processes. 
Really, typical particle masses are of order itlt ~ gT, therefore, for small g 
the Compton wave length Xcompt. ~ 1/^t, giving the particle size, and the 
Larmor radius, rL armor ~ l/(gH) l l 2 ~ l/g 2 T , determining the space range 
where the charged particle spectrum is formed, are both located inside a domain 
which is filled up by the field. Hence it follows that the field strengths of order 
(gH) 1 / 2 ~ g 2 T or stronger are of interest at high temperature and, in particular, 
the Sawidy mechanism gives rise such intense fields. 

Interesting mechanism to generate the hypercharge magnetic field due to the 
abelian anomaly was proposed in Ref. |7[|. At presen time it is not investigated 
in detail. So, in what follows we will not consider a consistent picture with this 
field. We will just assume that the strong external field Hy present when the 
transition had happened. 

1 In paper pp| the important cancellation of logarithmic terms entering the zero and the 
finite temperature parts of the EP was missed. This has been resulted in the incorrect value of 
the vacuum magnetic field spontaneously created at finite temperature. 
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Below, we shall study the influence of constant the usual magnetic []I8"|| , fllCi 
and the hypercharge magnetic 0, ||, 0, |l2j fields. Such the approximation is 
not artificial for strong fields. More definitely, the gradients of fields are negligible 
if the relation | VH/H \« H/m holds ||49|| . Here, m is a characteristic mass 
of the problem under consideration. It means that in strong fields the dominat 
effects are due to intensity of the fields. At high temperature, particle masses 
are of order m ~ gT. As it was pointed out in Ref. f| this approximation works 
well for the second order phase transitions and for the first stages of the first 
order ones when the bubbles are not large. On the other hand, in Refs. , f| it 
was argued that strong stochastic hypermagnetic fields are able to produce the 
baryon asymmetry of the universe. These problems are left beyond the scope of 
the present paper. 

2.2 Mechanisms of acting of the external fields on a vacuum 

The ways in which the magnetic and the hypermagnetic fields affect the vacuum 
scalar condensate are quite different. In the latter case, it is completely similar to 
the case of superconductivity, as it was investigated first in Ref. |L6|. In broken 
phase, the gauge field ( U(l) gauge field in the Higgs model and Z-boson field in 
the EW theory) is screened by its mass. This is the consequence of the interaction 
term ~ A^cf 2 presenting in the Lagrangian, is the scalar field. The influence of 
the external field is reduced to the increase of the vacuum energy and it manifests 
itself at tree level. In sufficiently strong fields the symmetry restoration happens 
and the gauge field mass Ma = g<p c becomes equal zero. For the critical field 
strength one has an estimate, H Y /2 ~ 0(M\), which shows that the restoration 
happens when the energy density of the external field equals to that of the scalar 
condensate ~ M\. This mechanism works both at zero and finite temperature. 
In the latter case, the critical value of Hy is decreasing when the temperature is 
increasing, as in superconductors. If the mass itlh < m z , the vacuum of the SM 
is the first type "superconductor" with respect to hypermagnetic field otherwise 



it is the second type one. This picture is determined at a classical level [|16 
The external hypermagnetic field delays the first order phase transition making 
it stronger, that is necessary for baryogenesis. This conclusion has been obtain in 



three approximation in Ref. ||. As it was discussed in detail in Ref. |Tj|, to have 
a standard baryogenesis scenario the ratio R = A<p c /T c of the order parameter 
jump to the critical temperature must be of order ~ 1.2 — 1.5 At zero external 
fields, one has R ~ 0.6 for run ~ 70 — 80 GeV. 

In the EW theory, the U(l) symmetry corresponding to electromagnetic field 
is retained in the broken phase. Therefore, there are no couplings of the scalar and 
the electromagnetic fields at tree level. Hence, one would expect that the external 
fields affect the vacuum condensate through the radiation corrections as it was 



discussed in Ref. [53]. However, in the nonabelian case the influence of fields 



on the scalar condensate is more complicate and the phase transition in strong 
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external magnetic fields can happen at tree level, because of the non-linearity 
of the field equations. The vacuum properties are determined by not only the 
scalar field condensate </> c but also the other order parameter - so-called P^-boson 

Actually, here one faces the 



condensate [pl| , [|43|| , |3ql (see also surveys |33 



situation with interacting order parameters. In the EW Lagrangian the next two 
terms, 



and ~ W+W~ , , 
field strength tensor, is the W-boson field 



enter, where F^ v is an electromagnetic 
The former interaction is due 



to a gyromagnetic ratio 7 = 2 inherent the nonabelian theories. It is crucial 
for exciting in strong constant mag netic fields H > H = M 2 /e in the W- 
boson spectrum of the tachyonic (unstable) mode p\ = p| + M 2 — eH, p 3 is a 
momentum along the field direction, which then is condensed owing to the self- 
interaction ~ (W+W~) 2 . The iy-boson condensate influences the scalar field 
at classical level (due to the second of the above written interaction terms). As 
a result, the scalar condensate is eliminated and the condensates of W- and Z- 



boson fields are formed p3[ , ||43|| , |pq| . The threshold of the phase transition is 
determined by the value of the Higgs boson mass. If m# < M w the vacuum 
is the first type "supercondustor" with respect to the magnetic field. In this 
case the homogeneous jy-condensate is formed [54|. For m# > M w the vacuum 
behaves as a second type "superconductor", therefore the lattice structure of the 
Abrikosov type formed by the W- and Z-boson fields is produced. This picture 
has been derived for arbitrary mass m# in Refs. |43|], |58| and for the special 
value m H = m z in Refs. |j3~6| , |32| . 



These vacuum properties have been determined in tree approximation when 
fermions do not contribute. But they do contribute in one-loop order, and, more- 
over, play a decisive role in the vacuum dynamics. In strong fields at high tem- 
peratures the influence of light fermions is increased, as the consequence of the 
term ~ H 2 log(T /rrif) entering the one-loop EP. Besides, other peculiarities must 
be taken into consideration. Namely, the EP contains additional T-dependent 
terms, such as ~ 4>lT 2 , ~ eHT(j) c , etc. (see more details in sects. 5, 7), which 
are generated in the one-loop and higher orders of perturbation theory. These 
terms make the picture of the field action involving. Therefore, to investigate 
symmetry behaviour numeric calculations must be carried out. 



2.3 Symmetry behaviour in magnetic fields 

As it is well known j|5| , j|6| , [[55] , at finite temperature side by side with the one- 
loop EP the correlation corrections described by the daisy diagrams should be 
added in the consistent calculations. Series of these diagrams are responsible for 
the long range effects and contain imaginary terms which cancel the imaginary 
part of the one-loop EP. As a result, the total EP is real at sufficiently high 
temperature. This important property is fulfilled also in the external fields when 
the contribution of daisy diagrams with unstable mode is allowed for [[[(J, |T . 



Now, we are going to describe symmetry behaviour at high temperature and 
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V(0 c ,T) 




Figure 1: Symmetry behaviour for the first order phase transition. T Cl and T c 
are upper and lower spinodal temperatures, T c is critiacal temperature. 



strong magnetic fields. In what follows we assume that only one type of the fields 
is applied. Remind that in the broken phase the component of Hy responsible for 
the magnetic field is unscreened and the form of the EP curve corresponds to both 
the magnetic and the hypemagnetic fields. In order to investigate the EW phase 
transition we shall consider the function V = V(h, B, 0) — V(h, B, 0) describing 
the symmetry restoration. Here, the dimensionless magnetic, h = H/H ,H = 
M^/e, and the scalar, (ft = <p c /6(0), fields and the inverse temperature, B = (3M W , 
are introduced, 5(0) is the EP minimum position at zero temperatures and fields 
(more details see in sects. 8, 9). The normalized order parameter <ft is changed 
from unit to zero that is convenient in numeric calculations. 

In fact, symmetry breaking (restoration) can be realized in two ways. By the 
first order phase transition with a nonzero jump of the order parameter A<p c ^ 
or by the second order one when the order parameter is changing smoothly. 

In Fig. 1 we depict symmetry behaviour at high temperatures for the first 
order phase transition. Notice that we consider the metastable electroweak min- 
imum which is separated from a global unbounded state disposed at large values 
of 2 by a wide and high potential barrier. The temperatures T C1 and T C2 are 
called the spinodal temperatures, which correspond the situations when bubles 
of the broken (T Cl ) and the restored (T C2 ) phases can not exist in the vacuum. 

In Fig. 2 we show symmetry behaviour for the second order phase transition. 
In this case baryogenesis can not be realized because of washing away the baryon- 



antibaryon asymmetry (see survey ||13||). 



As we mentioned before, at zero external fields the value of the parameter 
R ~ 0.6 and so baryogenesis is not possible in the SM. Our main task here is to 
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V(0 c ,T) 




Figure 2: Symmetry behaviour for the second order phase transition. The order 
parameter (p c (H,T) is increasing monotonically with temperature decreasing. 



investigate the situation with strong external fields been taken into consideration. 

The general picture of the phase transitions is independent of the specific 
external conditions considered. The external fields change the parameters (or the 
type) of the phase transitions. Besides, because of the different mechanisms of 
influence of the magnetic and hypermagnetic fields on the vacuum the conditions 
fixing the transition temperature are also different for these fields (for more details 
see sects. 8, 9). 

Now, let us describe the restored phase. For the magnetic field case, because 
of the presence of the tachyonic mode in the W-boson spectrum, the state = 
is unstable at zero temperature. To better understand the situation at T ^ 0, let 
us first discussed the one-loop case. To verify whether the vacuum is stable or 
not one has to consider the effective mass squared of the unstable mode defined 
as the sum of the ground state energy squared taken at zero momentum and 
the one-loop W-boson mass operator averaged over the ground state. It must be 
considered for the value of H condensed at high temperature. In particular, we 
have for the field condensed at the one-loop level (gH) l J 2 = (g 2 /2ir)T \$X§ : 
M 2 (H C , T) = H(H C , T, n = 0, a = +1) —gH c > 0. Thus, the vacuum stabilization 
is observed. The situation which takes place when the correlation corrections are 
included is discussed in sect. 11. 

In the restored phase the H^-bosons do not interact with the hypermagnetic 
field. So, no instabilities occur in this case. 

By investigating symmetry behaviour in either the hypermagnetic or the mag- 
netic field within the total EP inlcuding the contributions of all the SM particles 
we have determined by numeric computations that for all m# values considered 
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the increase in the field strength is resulted in the weaker (not stronger) first 
order phase transition. For H,Hy > 0.1 — 0.5 10 24 G it becomes of the second 
order at all. Thus, we come to the conclusion that baryogenesis does not survive 
at these external conditions. 

To better understand the role of fermions in symmetry behaviour let us adduce 
two terms of the asymptotic expansion of the one-loop EP in the limit of T — > 
oo, H — ► oo. The first one is the term ~ H 2 log^j. Due to this term the light 
fermions are dominant at high temperature. The second term can be derived 
from the expansion of the zero temperature part: ~ —eHm'j This term acts to 
make "heavier" the Higgs particles in the field. As a result, the second order 
temperature phase transition is stimulated due to strong fields. 

We would like to complete this section with the comparison of the described 
results with that of other approaches (for more details see sect. 8, 12). The EW 
phase transition in the hypermagnetic field has been investigated in one-loop ap- 
proximation to the EP in Ref. M and by the method combining perturbative 
results and lattice simulations in Ref. 0. These authors bacause of different 
reasons have skipped the fermion part of the EP and therefore had no possibil- 
ity to determine the form of the EP curve in strong fields at high temperature. 
In this respect our investigation filled up the gap existed. In our investigation 
the external fields have been taken into account exactly through Green's func- 
tions. Therefore, in particular, influence of fermions on symmetry behaviour was 
correctly reproduced. 



3. BOSON FIELD CONTRIBUTIONS TO V^{T, H,<f> c ) 

The Lagrangian of the boson sector of the Salam- Weinberg model is (see, for 
example, [BO]) 



L = --Ft, *r - - A G„ V G^ + (D^)+(D^) 



+ ^!($+$) _ ^($+$) 2 , (2) 



where the standard notations are introduced 



F^^d^At-d^ + ge^Al, 
D, = d, + \igAlr a + Ug'B,. 



The vacuum expectation value of the field $ is 



(3) 



< * >= f I ) • (4) 
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The fields corresponding to the W-, Z-bosons and photons, respectively, are 



1 



(5) 



To incorporate interaction with an external hypermagnetic field we add the 
term jHHy to the Lagrangian. The value of the macroscopic magnetic field 
generated inside the system will be determined by minimization of free energy. 
Interaction with classical electromagnetic field is introduced as usually by split- 
ting the potential in two parts: = A^ + A R , where A R describes a radiation 
field and A = (0,0, Hx 1 ,0) corresponds to the constant magnetic field directed 



along the third axis. We make use of the gauge- fixing conditions [33 



d,W^ ± ieApW^ =F i 9 ^-^ = C ± {x), (6) 

d^-^(g 2 + g y/ 2 M z = C z (x), (7) 

where e = gsin6 w ,tanO w = g' / g, (p^, (p z are the Goldstone fields, £, £' are the 
gauge fixing parameters, C ,C Z are arbitrary functions and <p c is a value of the 
scalar field condensate. In what follows, all calculations will be carried out in the 
general relativistic renormalizable gauge an d after that we set £, £' = 

choosing the unitary gauge. 

To compute the EP in the background magnetic field let us introduce 
the proper time, s-representation, for the Green functions 

oo 

G ab = -i J dsexp(-isG- iab ) (8) 



and make use the method of Ref . |59| , allowing in a natural way to incorporate 



the temperature into this formalism. A basic formula of Ref. |59| connecting 



the Matsubara Green functions with the Green functions at zero temperature is 
needed, 

+oo 

Gf(x, x'; T) = J2(-l) {n+[x])r7k Gt b (x ~ [x]Pu, x' - n(3u), (9) 

— oo 

where Gf 1 is the corresponding function at T = 0,(3 = l/T,u = (0,0,0, 1), the 
symbol [x] denotes an integer part of x±/ (3, — 1 in the case of physical fermions 
and cTfc = for boson and ghost fields. The Green functions in the right-hand side 



11 



of formula (|9[) are the matrix elements of the operators Gk computed in the states 
| x', a) at T = 0, and in the left-hand side the operators are averaged over the 
states with T^O. The corresponding functional spaces U° and U T are different 
but in the limit of T — > U T transforms into U°. 

The one-loop contribution into EP is given by the expression 

= ~Tr log G ab , (10) 

where G ab stands for the propagators of all the quantum fields W , (/> ,-■■ in the 
background magnetic field H. In the s-representation the calculation of the trace 
can be done in accordance with formula 

Tr log G ab = -J— trexp{-tsG^). (11) 
o 

Details of calculations based on the s-representation and the formula can be 
found in Refs. [[44]], [[0|]. The terms with n = in Eqs.((|), fllCf) give zero 



temperature expressions for the Green functions and the effective potential 
respectively. They are the only terms possessing divergences. To eliminate them 
and uniquely fix the potential we use the following renormalization conditions at 
H,T = ||: 

d 2 V((j),H) 1 



(12) 

H=0,<f>=8(0) = 0, (13) 

|h=o^=«5(o) = | m 2 |, (14) 



dH2 |H=O,0=<5(O) 2 



d 2 V{(f),H) 



where 1/(0, H) = V<® + + • • • is the expansion in a number of loops and (5(0) 
is the vacuum value of the scalar field determined in tree approximation. 

It is convenient for what follows to introduce the dimensionless quantities: 
h = H/H (H = M%/e),<f> = <j> c /5(0),K = m 2 H /M 2 , B = (3M Wj t = 1/B = 
T/M w , V = V/Hq and M w = §5(0). After reparametrization the scalar field 
potential is explicitly expressed through the ratio K, 

VW = y + irsin 2 ^(-0 2 + ^). (15) 

Notice that h in the case of the external hypermagnetic field is the component 
of hy which remains unscreened in the broken phase. In the restored phase, it 
will be convenient to work in terms of the initial fields and we will carry out the 
corresponding calculations later. 
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The renormalized one-loop EP is given by the sum of the functions 

V 1 = V (0) + V (1) (0, h, K) + cu (1) (0, h, K, r), (16) 
where is the one-loop EP at T = 0, which has been studied already in Ref. 



3l|. It has the form: 

V (1) =V«+V« (17) 

where 

v«= * [*v4 + |) + «'(-D + ^-^4 k ' 

- ift 2 ! og (2ft)] + ^|-2ft 2 + (ft 2 + h4?)log(h + 

+ (/i 2 - h4> 2 )log \h-4> 2 \}+ i^ah{<f) 2 - h)6(h - 2 ), (18) 



+ ^r l( /- 3 * + 2 )!o9|Ji 2- | -T^ (19) 

and wP) is the temperature dependent contribution to the EP given by the terms 
of formulae <g), (0) with n ^ 0. 

We outline the used calculation procedure considering the P^-boson contribu- 
tion as an example |T9| . 



OO 



4 1} = f / * e-^)[ 1+2C ° S2S ] £exp(MV/4 S ). (20) 

ZiTX J s sin 5 i 

1 

As Eq. (|T8|), this expression contains an imaginary part for h > <ft 2 appearing 
due to the tachyonic mode e 2 = p 2 + M 2 — eH in the VF-boson spectrum [Q . 
It can be explicitly calculated by means of the analytic continuation taking into 
account the shift s — > s — iO in the s-plane. Fixing (p 2 /h > 1 one can rotate 
clockwise the integration contour in the s-plane and direct it along the negative 
imaginary axis. Then, using the identity 

1 oo 



oo 



sinhs p=0 



and integrating over s in accordance with the standard formula 

b , ,b. 

dss' 1 x exp(- 

n 
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h h 

f dss" 1 - 1 exp( as) = 2(-) n/2 K n (2Vab), (22) 

J s a 



a, b > 0, one can represent the expression (p0[) in the form 



Rew$ = -4-4(3wo +wi - w 2 ), (23) 

7T £5 



where 



oo oo 

= E E -^iW; a; p = (0 2 + /i + 2P/1) 1 / 2 (24) 

p=0n=l n 

00 7/ 

Wl = E-^i(^),2/=(0 2 -^) 1/2 (25) 



and in the range of parameters (jr < h after analytic continuation 



^-jE-^MIl/l), (26) 

Z n=l n 



oo 

^ 2 = E -K,{nBz),z = (0 2 + /i) 1 / 2 , (27) 

i<r n (x), l^i (x) are the MacDonald and Bessel functions, respectively. The imagi- 
nary part of uj^ is given by the expression 

h 00 I v I 

Imwt = -2a- E — MnB \ y |), (28) 

Ji(x) is Bessel function. As it is well known, the imaginary term of the EP is 
signaling the instability of a system. In what follows, we shall consider mainly 
symmetry behaviour described by the real part of the EP. As the imaginary part 
is concerned, it will be cancelled in a consistent calculation including the one-loop 
and daisy diagram contributions to the EP. 



Carrying out similar calculations for the Z- and Higgs bosons, we obtain : 



= - 6 - £ J *m K *^ ( 29 ) 

7r cos z u w n z B z cos v 

Re ^ = iaz JS,y,(«B 1 1 1) }■ (30) 

n=l 

where the variable t = {K w ( 3 ^ 2 ~ 1 - > )] 1 ^ 2 at 3</> 2 > 1 and the series with the function 
Y2(x) has to be calculated at 3(f) 2 < 1. The corresponding imaginary terms are 
also cancelled as it will be shown below. 



The above expressions QT7D, (|23|), ( |30D will be used in numerical studying 
of symmetry behaviour at different H, T. Notice the cancellation of a number of 
terms entering the zero-temperature part given Eqs. fllTD and T-depended one. 
This fact has a general character and was used in checking of the correctness of 
calculations. 
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4. Fermion contributions to VM(H,T, 



The fermion one-loop EP in magnetic fields is well studied [53], [29], [63]. To find 
this explicit form at finite temperature let us consider the standard unrenormal- 
ized expression written in the s-representation 

= JL £ (-1)" / ^e-K^ 2 ™ 2 /4,) (ejffs)co ^ (ejffs); (31) 

m/ is the fermion mass. Here, we have incorporated the equation @ to include a 
temperature dependence. In what follows, we shall allow for the contributions of 
all fermions - leptons and quarks - with their present masses. Usually, considering 
symmetry behaviour without field one restricts himself by a t-quark contribution, 
only. But in the case of external fields applied this is not a good idea, since the 
dependence of on H is a complicate function of the ratio rrij/eH. So, at 
some fixed values of H, T the different kind dependencies on H will contribute for 
fermions with different masses. Hence, a very complicate dependence on the field 
takes place in general. We shall include this in a total by carrying out numerical 
computetions and summing up over all the fermions. Now, separating a zero 

+0O oo 

temperature contribution by means of the relation — 1 + 2 J2 and introducing 

— oo 1 

the parameter Kf = m'j/M^ = G\ ukaw j ' g 2 ', we obtain for the zero temperature 
fermion contribution to the dimensionless EP, 



/ 



a v^/ 2 | q f | h 

- ^p^iogr l( ^) + (2C(-i)-i)^ 



+ \ K ^ 1 +(\Ky-\K l \ qt \h4?)\o i 2 -^-], (32) 



3 3 

where qt is a fermion electric charge, the sum J2 = J2 (leptons) + 3 J2 (quarks) 

f /=i f=i 
counts the contributions of leptons and quarks with their electric charges. The 



Ti function is defined as follows (see Refs. |29[, [p3|]): 



J 1 1 

logTi(x) = J dy\ogT(y) + -x(x - 1) - -slog(27r). (33) 



15 



w \ 

/ \ 



^' w >v 

/ \ 



/ 



/ \ / 

I 1 + 1 

\ / \ 

X / \ y 



\ I 

I + I 

/ • 



Figure 3: The Higgs field daisy diagrams giving contribution to the effective 
potential. Blobs stand for the neutral scalar field polarization operator calculated 
at zero momentum. 



The finite temperature part can be calculated in a way described in the pre- 
vious section. In the dimensionless variables it looks as follows: 



UJr 

7T 

" / p=0n=l 



+ «2P + 2)A + K,^" + + K ^ Bn)] X (34) 

Again, a number of terms in Eqs. ( p2| ) and ( p4| ) are cancelled in the total, as in 
the bosonic sector. 

These two expressions will also be used in numeric investigations of symmetry 
behaviour. 

5. CONTRIBUTION OF DAISY DIAGRAMS 



It was shown by Carrington |Jd) that at T ^ the consistent calculation of the 
EP based on the generalized propagators, which include the polarization operator 
insertions, requires that daisy (ring) diagrams have to be added simultaneously 
with the one-loop contributions. These diagrams essentially affect the phase 
transition at high temperature and zero field P5(| , ||46|| , f55|j . Their importance at 



T and H ^ was also pointed out in Refs. [10], [11 



As it is known |Tj|, ||45|| , the sum of daisy diagrams describes a dominant 
contribution of long distances. It is important when massless states appear in 
a system. So, this type of diagrams has to be allowed for when a symmetry 
restoration is investigated. To find the correction V ring (H, T) at high temperature 
and magnetic field the polarization operators of the Higgs particle, photon and 
Z-boson at the considered background should be calculated. Then, V ring (H,T) 
is given by a series depicted in figures 3, 4. Here, a dashed line describes 
the Higgs particles, the wavy lines represent photons and Z-bosons, the blobs 
correspond to the polarization operators in the limit of zero momenta. As also it 
is known |45[| , 4TJ, in order to calculate the contribution of daisy diagrams the 
limiting expressions of the polarization operators (k, T, H) at zero momenta, 
noo(^ = 0, T, H), are to be substituted. This limit, called the Debye mass, can be 
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• o • • 

Figure 4: The photon and Z-boson daisy diagrams giving contribution to the 
effective potential. Blobs stand for the polarization operators of the fields calcu- 
lated at zero momenta. 




computed from the EP of the special type. The latter fact considerably simplifies 
our task. 

Now, let us turn to calculations of V r i ng (H,T). It is given by the standard 
expression fll5f , |45[| , [[46] , |i~7f : 



V ring = -7^>{[M 2 (0) + n 00 (0)] 3/2 - Af 3 (0)}, (35) 

where trace means the summation over all the contributing states, M(0) is 
the tree mass of the corresponding state. The functions n 00 (0) are: n 00 (0) 
=U(k = 0,T,H) for the Higgs particle; n 00 (0) = n 00 (£: = 0,T,H) - the zero- 
zero components of the polarization functions of gauge fields in the magnetic 
field taken at zero momenta. The above contributions are of order 
in the coupling constants whereas the two-loop terms have order ~ g A (X 2 ). For 
Hoo(0) the high temperature limits of polarization functions have to be substi- 
tuted which are of order ~ T 2 for leading terms and ~ g<fi c T, (gH) l l 2 T((f) c /T « 
1, {gH) l ' 2 /T « 1) for subleading ones. 



5.1 The polarization function of scalar field 

For the next step of calculation, we remind that the effective potential is the 
generating functional of the one-particle irreducible Green functions at zero ex- 
ternal momenta. So, to have 11(0) we can just calculate the second derivative 
with respect to (j) of the potential V^'(H, T, (ft) in the limit of high temperature, 
T » <f),T » (eH) 1 / 2 , and then set = 0. This limit can be calculated by 
means of the Mellin transformation technique (see, for instance, |]19[) and the 
result looks as follows: 

+ if (3*2 - m)T> - - Vr - f( 3A ^:* 2(0) pr] 

6 cos 6 6 6 2 
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+ L eHT {Ul-eHY'\ (36) 
271 4 

3 3 

The parameter Cj = J2 G\ + 3 J2 Gf q determines the fermion contribution of 

i=l i=l 

order ~ T 2 having relevance to our problem. We also have omitted ~ T 4 contri- 
butions to the EP. The terms of the type ~ log[T '/ f((j), H)] cancel the logarithmic 
terms in the temperature independent parts fllTf), (|31"|). Considering the high tem- 
perature limit we restrict ourselves to linear and quadratic in T terms, only. 

One else important expression, which also should be taken into account, is 
the linear in H term of the zero temperature EP which looks as follows: 

vff(H,4> c )/H 2 = -|/X>/ I IfH I • (37) 

It significantly affects symmetry behaviour and contributes to the Debye mass in 
strong fields. 

Now, differentiating these expressions twice with respect to <fi and setting 
= 0, we obtain 



n*(o) 



d 2 V^((f) } H,T) 

al2 U=° 



, 6e 2 3e 2 
■ 6A H ry—— + 

V sin 2 2/9. . 



24/? 2 V sin 2 29 w sin 2 9 V 

-Y.[- w -\i,»\k. 



+ J eH t'\ e 2 (3^«-lb- (38) 
8ir sin 9 w p 2 2 

The terms ~ T 2 give standard contributions to temperature mass squared com- 
ing from the boson and fermion sectors. The if-dependent term is negative since 
the difference in the brackets is 3v / 2C( — 1> \) ~ ^ — — 0? 39. Formally, this may 
result in the negativeness of 11(0)0 for very strong fields (e#) 1/2 > T. But this 
happens in the range of parameters where asymptotic axpansion is not applica- 
ble. Substituting expression (|38D into Eq. (|35|) we obtain (in the dimensionless 
variables) 

VL, = ~{(^K + n,(0)} 3/2 + ^K^f*. (39) 

As one can see, the last term of this expression cancels the fourth term in Eq. 
(|36|). which becomes imaginary at 3(f) 2 < 1. This is the important cancellation 
preventing the infrared instability at high temperature. 

Notice that Eq. ([3(]) contains other term (the last one) which becomes imag- 
inary for strong magnetic fields or small 2 . It reflects the known instability in 
the W-bosoia spectrum ppf , |3"T[ . 
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5.2 The Debye masses of photons and Z-bosons 

To find the if -dependent Debye masses of photons and Z-bosons the following 
procedure will be used. We calculate the one-loop contributions to the EP due to 
the P^-bosons and the fermions in a magnetic field and some "chemical potential" , 
/i, which plays the role of an auxiliary parameter. Then, by differentiating them 
twice with respect to fi and setting fi = the mass squared m 2 D will be obtained. 
Let us outline that in more detail for the case of fermion contributions where the 
result is well known. 

The temperature dependent part of the one-loop EP of constant magnetic 
field at a non-zero chemical potential induced by an electron-positron vacuum 



polarization is [63]: 



1 00 yjcj — /9 2 / 2 
^eL = ^E(-l) m y m 2 s)(eHs)coth(eHs)cosh((3l^ (40) 

1=1 o 

where m is the electron mass, e = gsin9 w is the electric charge and s-representation 
is used. Its second derivative with respect to fi taken at fi = can be written in 
the form 



= ^{3 2 ^ £(-l) i+1 / -exp(-m\s - (3 2 l 2 /As)coth(eHs). (41) 

Expanding coth(eHs) in series and integrating over s in accordance with formula 
(|22|) we obtain in the limit of T » m,T » (e#) 1/2 : 

E("l) m &i(^) + l^^-K^l) + ■••]. (42) 
pi dm 

Series in / can easily be calculated by means of the Mellin transformation (see 
Refs. |0, jn3|, [0). To have the Debye mass squared it is necessary to dif- 



ferentiate Eq. (fllD with respect to (3 2 and take into account the relation of the 
parameter [i with the zero component of the electromagnetic potential : /i — > icAq 



17]. In this way we obtain finally, 



ml = g 2 sin 2 9 w (^- - -^m 2 + 0((m/5) 2 , (eH{3 2 ))). (43) 



This coinsides with the known result calculated from the photon polarization 
operator [|62]| . As one can see, the dependence on H appears in the order ~ T~ 2 . 



To find the total fermion contribution to m D one has to sum up the expression 



( f43D over all fermions and substitute their electric charges. 

To calculate m 2 D for Z-bosons it is sufficient to account of the fermion coupling 
to Z-field. It can be done by substituting jj, — > i(g/2cos9 w + gsin 2 9 w ). The result 
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differs from Eq. fl43| ) by the coefficient at the brackets in the right-hand side 
which has to be replaced, g 2 sin 2 9 w — > ^ 2 ( 4eo ] 2e + tang 2 6 w ). One also has to add 
the terms coming due to the neutral currents and the part of fermion-Z-boson 
interaction which is not reproduced by the above substitution: 

+ (44) 

8cos z m 



Now, let us apply the above procedure to obtain the VT-boson contribution. 
The corresponding EP (temperature dependent part) calculated at non-zero T, /z 
looks as follows, 

V£> = -f^ £ ]%xp(-m 2 s - 1 2 P 2 /4s)[ , , f - . + Asinh(eHs)]cosh(M. 
o7r J s z sinh{eHs) 

(45) 

All the notations are obvious. The first term in the squared brackets gives the 
triple contribution of the charged scalar field and the second one is due to the 
interaction with a VT-boson magnetic moment. Again, after differentiation twice 
with respect to /i and setting \i = it can be written as 

d 2 V^ eR d ™7ds . m 2 s l 2 (3 2 eH.. 3 A . , . 

- Z = — j/r^V / —exp( — - A ) —-+Asinh(s)}. (46) 

<V 2vr 2M a/3 2 ^i s FK eH As n sinh{s) v n K ! 

Expanding sinh~ 1 s in series over Bernoulli's polynomials, 

i — s oo r> 

— Ett(-2^, (47) 



sinks s £^f Q k\ 



and carrying out all the calculations described above, we obtain for the W- boson 
contribution to m 2 D of the electromagnetic field 



m 

0(m 2 /T 2 ,(^r^tf/T 2 ) 2 )]. (4* 



m z? = 3g 2 sin 2 w [-T 2 - —T{m 2 + gsin6 w H) 1/2 - -^r(gsin8 w H) 

6 Zll 07T Z 



Hence it follows that spin does not affect the Debye mass in leading order. Other 
interesting feature is that the next-to-leading terms are negative. 

The contribution of the H^-boson sector to the Z-boson mass m 2 D is given 
by expression (f||) with the replacement g 2 sin 2 6 w — ► g 2 cos 2 6 w . Summing up the 
expressions (ffi^) and (HSj) and substituting them in Eq. fl3"5D, we obtain the photon 
part V^] n9 , where it is necessary to express masses in terms of the vacuum value 
of the scalar condensate <p c . In the same way the daisy diagrams of Z-bosons 
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V r z ing can be calculated. The only difference is the mass term of Z-field and the 
additional term in the Debye mass due to the neutral current ~ v^^vZ^. These 
three fields - <p, 7, Z- which become massless in the restored phase, contribute 
into V r i ng (H, T) in the presence of the magnetic field. At zero field, there are also 
terms due to the jy-boson loops in Figs. 3, 4. But when H 7^ the charged 
particles acquire masses ~ eH and these daisies can be neglected. 



5. 3 Daisy diagrams of the tachyonic mode 

A separate consideration should be spared to the tachyonic (unstable) mode in 
the IV-boson spectrum: p\ = p\ + M 2 — eH. First of all, we notice that this mode 
is excited due to a spin interaction and it does not influence the Goo(&) component 
of the W-bosoia propagator. Secondly, in the fields eH ~ M 2 , the mode becomes 
the long range state. Therefore, it should be included in V r i ng (H,T) side by 
side with other considered neutral fields. But in this case it is impossible to 
take advantage of formula (|35|) and we return to the initial EP containing the 
generalized propagators. 

For our purpose it will be convenient to make use of the generalized EP written 
as the sum over the modes in the external magnetic field [[17]], [18| : 

K$ = ^ E / £ E + <U P3 + n(T, H))), (49) 



/= — oo_„ ra=0,cr=0,±l 



where LU t = ?f , el = p\ + M 2 + (2n + 1 - 2a)eH and U(H, T) is the radiation 
mass squared of W^-bosons in magnetic field at finite temperature. Denoting as 
Dq l {pz, H.T) the sum uf + e 2 , one can rewrite eq. (f49l) as follows: 



+ SE + [^{log[l + (ujf+pt + M 2 w -eH)^U(H,T)] 

+ E log[l + D (elH,T)U(H,T)]}. (50) 

Here, the first term is just the one-loop contribution of H^-bosons, the second one 
gives the sum of daisy diagrams of the unstable mode (as it can easily be checked 
by expanding the logarithm in a series). The last term describes the sum of the 
short range modes in the magnetic field and should be omitted. 

Thus, to find V^f g table one has to calculate the second term in Eq. (|50"D. In 
the high temperature limit we obtain: 

P 11 

V ZT hle = ~ eH + n(^T)) 1 / 2 - {Ml - eH)^}. (51) 
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By summing up the one-loop EP and all the terms V r i ng , we arrive at the total 
consistent in leading order EP. 

Let us mention the most important features of the above expression. It is seen 



that the last term in Eq. (51) exactly cancels the "dangerous" term in Eq, 
So, the EP is real and no instabilities appear at sufficently high temperatures 
when Il(iJ, T) > M 2 — eH. To make a quantitative estimate of the range of 
validity of the total EP it is necessary to calculate the P^-boson mass operator 
in a magnetic field at finite temperature and hence to find U(H,T). This is a 
separate and enough cogent problem which is considered in detail in Ref. [g6 . 
Here, we only adduce the result of IT(if, T) calculations: 

n unstMe (H,T) = <n = 0,a = l\U c fl h u ar9ed \n = 0,a = l> 

= a[12,33(e#) 1/2 T + *4(e#) 1/2 T], (52) 

where the average value of the mass operator in the ground state of the W- 
boson spectrum | n = 0, o = +1 > was calculated. This expression has been 
obtained in the limit eH/T 2 « 1,B — M w /T « 1. Side by side with the 
real part responsible for the radiation mass squared the expression ( |52|) contains 
the imaginary one describing the decay of the state. The latter term is small as 
the former one is compared and of order of the usual damping constants at high 
temperature. So, ImU(H, T) can be ignored in our problem. The radiation mass 
squared is positive. It acts to stabilize the spectrum. At H = no screening is 
produced in one-loop order, as it should be at finite temperatures for transversal 
modes |TjJ. Thus, we come to conclusion that at sufficiently high temperature 



the effective iy-boson mass squared = M 2 — eH + H(H, T) is positive and 



no conditions for iy-boson condensation discussed in Refs. jJ3|, are realized. 


6. RESTORED PHASE IN THE EXTERNAL MAGNETIC 
FIELDS 

Having calculated the EP at 7^ we are able to determine the kind of the EW 
phase transition for different m#, h. That will be done in the next sections. Here, 
we are going to describe in more detail the restored phase with the hypermagnetic 
and magnetic fields. Let us consider first the former case. 

To describe more precise the restored phase one has to calculate radiation 
corrections to the external field Hy at high temperature. Before doing that let 



2 Express (52) disagrees with the corresponding one of Ref. [Q where the average value of 
the gluon polarization operator in an abelian chromomagnetic field was calculated in weak field 
approximation and 11(77, T) has been found to be zero. Most probably, the discrepancy is the 
concequence of the calculation procedure adopted by these authors when the gluon polarization 
operator was calculated at zero external field and then its average value has been calculated in 
the state | n = 0, a = +1 > . Our expression is the high temperature limit of the mass operator 
which takes into account the external field exactly. 
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us remind that at = this field is completely unscreened whereas the non- 
abelian constituents of the electromagnetic and Z-fields are screened at scales 
I > (<? 2 ^) 1 by the magnetic mass. Remind also that we are investigating the 
separate influence of the external fields. This means that in the covariant deriva- 
tive describing interaction with the external field Hy in the restored phase one 
should maintain the U(1) Y term only: = + \gB e * 1 . We set the potential 
as before, B* xt = (0, 0, H Y x\ 0). 

In the restored phase P^-bosons do not interact with Hy- The field dependent 
part of the EP V(<ft = 0, Hy, T) is non-zero due to the contributions of fermions 
and scalars. However, the fermion part depends logarithmically on temperature 
(~ ( g / 2 ) HylogT /T ) and can be neglected as compared to the tree level term 
\H Y . This is not the case for the scalar field whose contribution to the one-loop 
EP is 



v^(h y ,t)= - ( g, W 2R y 

sc ^ £1Y ' 1 ) 247T 2 



ln(T/T ) 



+ (W/2)H Y )^T 

07T 

The term logarithmically dependent on T can again be neglected but the linear 
in T part should be retained. Since "hyperphotons" are massless in the restored 
phase we also include the contribution of the corresponding daisy diagrams: 

V rest 9 ored( H Y,T) = ~ [| ( 9 > /2) 2 T 2 + m% - MBM^L - _L {g > /2 )Hy] ^ , 

" (54) 
where m 2 D = ^g' 2 T 2 Yf is the sum over the fermion contributions to the 

Debye mass of the "hyperphotons" , Yf are the hypercharges of R— and L— lep- 
tons and quarks. Both these expressions have been calculated in a way described 
in previous sections. 

For convenience of numerical investigations we express Eqs. fl53p and (1531) in 
terms of the dimensionless variables h, B: V(H Y , T) restored = (H ) 2 v restor ed(h, B), 

1 h 2 a h*' 2 

Vrestore d (h,B) - 3 ^ + £ 

1 a r 7 Ana h 1 / 2 h . . 

_ 1 !■>■- (55) 



3B l 6cos 2 6B 2 2y/2nBcosO 167i 2 cos 2 9 i 



where a = e 2 /47r and hy = h/cos6. 

Now, let us turn to the magnetic field case. In accordance with our approach 
we set H Y = 0. It will be important for what follows to remember recent results on 
obeservation of the gluon magnetic mass in lattice simulations that was found to 
be of order m mag ~ g 2 T (as it has been expected from nonperturbative calculatios 
in quantum field theory jTj|, [T3J]). The mass screens the nonabelian component 
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of the magnetic field at distancies / > I 

m ^ (s^-^ 1 ) ^ but insid.6 th.6 spcicc domain 
may exist and affect all the processes at high temperatures, as it was 
discussed in sect. 2.1. This fact has not been taken into consideration in a number 
of investigations of the EW phase transition. In particular, in Ref. |34j (as in Ref. 
[pT| ) the field strength generated at finite temperature was erroneously estimated 
as coinsiding with that at zero temperature. 

The magnetic field in the restored phase may be homogeneous or not depen- 
dently on the stability of the charged boson spectrum in the field calculated with 
radiation corrections included. That can be checked in accordance with formulae 
fl5T|) and (|52|) at M w = 0. If the effective mass squared M^ff (H, T) is positive, the 
perturbative vacuum is stable and the external field is homogeneous otherwise it 
is unstable and a lattice structure has to be generated due to the evolution of the 
instability. We shall see below, for the first order phase transition (which is the 
main topic in the present paper) M^j (H, T, C ) is positive in the minimum of the 
EP when the symmetry restoration happens. Therefore, we will not investigate 
in detail the structure of the restored phase for different H restricting ourselves 
by the case of the constant field, as in the broken phase. 

To make a link between studies of symmetry behaviour in the external hy- 
permagnetic field and the previous results for the usual magnetic field p4f , |L7| 



we notice that in the broken phase Hy and H are connected by the relation 
H = Hy cos 6. So, all investigations, dealing with symmetry behaviour in a mag- 
netic field at high temperature, are relevant in the case of Hy in the respect of 
the form of the EP curve at different T, Hy. The hyper charge field influences the 
scalar field condensate at tree level and acts to restore symmetry. That was the 
reason why it has been taken into account in the lower order in Refs. ||, H. 
But, as it will be shown below, for strong fields and heavy itlh the form of the 
EP curve in the broken phase is very sensitive to the change of the parameters. 
Moreover, it is strongly depended on the correlation correction contributions. So, 
to have an adequate picture of the EW phase transition symmetry behaviour with 
the daisy diagrams included has to be investigated. 



7. HIGH TEMPERATURE EXPANSION OF THE EF- 
FECTIVE POTENTIAL 

The most of investigations dealing with symmetry behaviour at high temperature 
used the limiting form of the EP at T — > oo. It will be of interest to compare the 
results obtained in two cases - for the complete EP and for the asymptotic one. 

First, let us adduce the high temperature limit of the sum of terms describ- 
ing the contributions of the scalar field at zero and finite temperature, the H- 
independent contribution of jy-bosons, Z-bosons as well as the contributions of 
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the second terms of corresponding daisy diagrams, Eq. (|39|): 

W = 8^^ (0 a + ^(M-B-) 2 -^)) (56) 

+ l2V^ (302 - 1)2 ^W + 602 - C(302 ~ 1)2] (5?) 
a,^^ n , tt 2 2 27T0 3 1 



7r L 12S 2 v r / 2cos 2 6B 2 3 cos 3 #5 J ' 

C= 0.5772 is Euler's constant. The first terms of daisies dependent on the Debye 
masses should be added separately. 

The high temperature asymptotic of the fermion sector u)f looks as folows: 

°* = ~D-yf + >Hs^)-^) (58) 

The term 47ra in Eq. (|58]) is introduced in order to correct the infrared infinity 
in the field. It effectively accounts of the fermion temperature mass. Just due to 
this infrared singularity the light fermions dominate in strong fields. 

Now, let us wright down the high temperature limits of the TV-boson sector 
(expressions ujq, lo\ — cj 2 ) Eq. (p3|) . Instead Re uj^ 1 one has to substitute the 

SUm U spin + Uq = J w . 

Uspin = -^[(0 2 + h)log(<p 2 + h)-(<p 2 -h)log\<p 2 -h\ (59) 
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Uo = J i [t k {l09 ( h + %i* B ~ c) - ( h + ^ V - c) m 



KW ^(0 2 + hf' 2 - —(h + 2 ) 1/2 + —h 2 - U 2 h - -< 



3B 2 3B XT ' B y T ' 16 8 16 
Here, the term due to the unstable mode is omitted since it is cancelled by the 
second term of daisy diagrams generating by the unstable mode. As it is seen, a 
lot of terms from expresion V z ^ at zero temperature are cancelled in the total. 
The term with log(Bft) is 0-independent and therefore inessential when symmetry 
behaviour is investigated, ft marks the normalization point. 

Now, for completeness, let us write down explicitely the contributions of the 
7- and Z- daisy diagrams: 

V r l ng = -^[n,(h,B)f 2 , (61) 
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where 

As the contribution of the Z boson daisy diagram we get (taking into account 
the subsitution e 2 — > 5 f2 ( 4eo 1 s 2 6) + tan 2 6) and addind other necessary terms): 

= -^[n,(A,S) + -^-] 3 / 2 + (63) 
y 3i5 cos^fl SB cos 6 U 



where 



8cos 2 9 [ 6B 2 2nB ) 



SB^cos^t) 



in 2 



The term, containing y K comes from the one-loop diagram of the Higgs field. 
Remind, for the case of the asymptotic EP the last term of Eq. (^3) has already 
been included in Eq. (|56|). 

8. SYMMETRY BEHAVIOUR IN STRONG HYPER- 
MAGNETIC FIELD 

Let us investigate the EW phase transition in the hypermagnetic field for different 
values of m# . It can be done by considering the Gibbs free energy in the broken, 
G br oken{H ext ,(j) c ,T), and the restored, G restored (H^ xt ,T), phases @, [§: 

G broken = V{<P,H,T)-HH e Y xt cose, (65) 

Grestored = V (0, Hy ,T) — HyHy 1 ■ 

By minimization of these equations the fields H and Hy generated in the vacuum 
have to be expressed throught Hy*. The first order phase transition can be 
determined within two equations: 

G r estored{HY Xt ,T,0) = Gb ro ken{H eXt ,T, <j)(H eXt ) c ) , (66) 
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describing the advantage of the broken phase creation, where 4>{H) C is a scalar 
field vacuum expectation value at given H, T, which has to be found as the 
minimum position of the total EP, 

dv(H,T^ c y° tal = o 

d(p c 

Hence the critical field strength can be calculated. In this expression and below 
we use for brevity H instead H ext . 

Having obtained the EP in the restored phase, the one-loop EP described by 
formulae (|lj]), (|23|) , (p9|)-(^4|) and the daisy diagram contributions V r i ng we are 



going to investigate symmetry behaviour. We shall present the results in two 
stages. First, we shall consider the total EP as the function of 2 at various fixed 
H, T, K and determine the form of the EP curves in the broken phase. In this 
way it will be possible to select the range of the parameters when the first order 
phase transition is realized. After that the temperature T c at given field strength 
Hy will be found. 

As usually ||33|| , to investigate symmetry behaviour we consider the difference 
V = Re[V(h, 0, K, B) - V(h, = 0, B)} which gives information about the sym- 
metry restoration. We will present the results for two cases: 1) for the total EP 
in the broken phase; 2) for the high temperature limits of it. Then we will make 
a comparison. 

In what follows, it will be also convenient to express the conditions of the phase 
transition in terms of the dimensionless variables h,B,<p, taking into account the 
relation hy = h/cos8. Then, the Gibbs free energy 

G broken (h ex \ 0, B) = y + v'(h, 0, B) - hh ex \ (68) 
has to be expressed in terms of h ext through the equation 

h»* = h+ dv ' {h >*> B \ (69) 
oh 

where v describes the one-loop and daisy diagram contributions to the EP. The 
phase transiton happens when the condition 

h 2 

—tan 2 6 = v restored (h, B c ) - v' broken {h, C , B c ) (70) 

holds. The function v' restored is given by Eq. (|5o]). We also have substituted the 
field h ext by h. 

The results on the phase transition determined by numeric investigation of 
the total EP are summarized in Table 1. 
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h 


K 


T c (GeV) 


<pc{h,T c ) 


<P 2 c(h,T c ) 


R 


M w {h,T c ) 


n m 

U.Ul 


U.oO 


1 Dfi A 7 


U.OU100Z 


u.uy i 


u.oyoyy 


U.OZ 1 ZOO 


0.01 


1.25 


122.21 


181659 


0.033 


36567 


230086 


0.01 


2 


145.56 


0.094868 


0.009 


0.16033 


0.186168 


0.1 


0.85 


108.58 


0.275681 


0.076 


0.62459 


0.245186 


0.1 


1.25 


123.54 


0.130384 


0.017 


0.25963 


0.112721 


0.1 


2 


148.39 


0.031623 


0.001 


0.05242 


0.126315 


0.5 


0.85 


108.89 


0.248998 


0.062 


0.56253 


0.49938 i 


0.5 


1.25 


second 


order 


phase 


transition 




0.5 


2 


second 


order 


phase 


transition 





Table 1. 

In the first column we show the hypermagnetic field strength in the broken 
phase (in dimensionless units). In the second and third ones the mass param- 
eter K = m 2 H /M 2 and the critical temperature of the first order phase transi- 
tion are adduced. Next two columns give the local minimum position <f) c (H, T c ) 
and its squared value at the transition temperature. The last two columns fix 
the ratio R = 246GeV(j) c (h, T c )/T c , determining the advantage of baryogenesis, 
and the P^-boson effective mass calculated in the local minimum of the EP at 
the corresponding field strength and the transition temperature. The parame- 
ter M w (h } T c ) = [(§0 C (>,T C )) 2 - eH + U(H, T c , (j) c )] 1 / 2 /M w is the dimensionless 
W-boson mass with the radiation corrections included. 

As it is seen, the increase in h makes the phase transition weaker (not stronger 
as it was expected in Refs. , || by analogy to superconductivity in the external 
magnetic field). The ratio R is less then unit for all the field strengths, wherease 



the baryogenesis condition is R> 1.2 — 1.5 ||13|| . Thus, we come to the conclusion 
that external hypermagnetic field does not make the EW phase transition strong 
enough to produce baryogenesis. Moreover, for strong fields the phase transition 
is of second order for all the values of K considered. 

Let us continue the analysis of data in the Table 1. For the field strengths 
H > 0.1— 0.5Hq(Hq = M 2 /e) the phase transition is of second or weak first-order. 
The W-boson effective mass squared (in dimensionless units) M 2 ((p c ,h,B c ) = 
(pl(h, B c ) — h + U(h, B c ) is positive for h = 0.01 and h = 0.1. Therefore, the local 
minimum is the stable state at the first order phase transition. For stronger fields, 
when the second order phase transition happens, the effective iy-boson mass 
becomes imaginary. This reflects the known instability in the external magnetic 
field which exhibits itself even when the radiation mass of the tachyonic mode is 
included. But it does not matter for the problem of searching for the strong first 
order phase transition in the external field investigated in the present paper. The 
instability has to result in the condensation of W- and Z-boson fields at high 
temperature. However, this does not change the type of the phase transition. 
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In Refs. 0, ||, ||, |HJ it has been concluded that the strong hypermagnetic 
field increases the strength of the first order phase transition and in this case 
baryogenesis survives in the SM. Our results are in obvious contradiction with 
this conslusion. To explane the origin of the discrepancies let us first consider 
Refs. [[J, [§] where a perturbative method of computations has been applied. 
These authors, in studying of the EW phase transition, have allowed for the 
influence of the external field at tree level, only. That corresponds to the usual 
case of superconductors in the external magnetic field, and, as a consequence, they 
observed the strong first order phase transition. In fact, the type of the phase 
transition was just assumed, since no investigations of the EP curve with all the 
particles included for different Hy, T have been carried out. In the former paper 
there was the qualitative estimate of the field effect, whereas in the latter one the 
quantitative analysis in one-loop approximation for the temperature dependent 
part of the EP has been done. Actually, in both these investigations the influence 
of the external field was reduced to the consideration of the condition ( [70|) fixing 
the transition temperature in the hypermagnetic field. The role of fermions and 
H^-bosons in the field was not investigated at all. However, as we have seen, the 
fermions (heavy and light) are of paramount importance in the phase transition 
dynamics. Just due to them the EW phase transition becomes of second order 
in strong fields (for the values of K when it is of first order in weak fields). 

In Refs. ||, |67| the phase transition was investigated by the method com- 
bining the perturbation theory and the lattice simulations. As the first step in 
this approach the static modes only are maintained in the high temperature La- 
grangian. The fermions are nonstatic modes and decoupled. The only fermion 
remainder is the t-quark mass entering the effective universal theory |T3|], M, 



57] . So, no fermion features in the external fields and hence no information 
about the form of the EP curve could be derived in this way. In our analysis, it 
has been observed that not only heavy but also light fermions are important in 
strong external fields, as one, in particular, can see from the term H 2 logT /uif 
of the one-loop EP. At high temperature it influences symmetry behaviour con- 
siderably. Actually, for various field strengths the fermoins with different masses 
are dominant and we have allowed for all of them. Besides, we have taken into 
consideration all the daisy correlation corrections in the external field that also 
affects symmetry behaviour. 

We would like to stress that our perturbative results for the values of K ~ 
0.8 — 0.9 are reliable. They are in agreement with nonperturbative analysis at 
zero field. The external field is taken into consideration exactly. For these mass 
ran ~ 75 — 80 GeV we observed the change of the first order phase transition into 
the second order one with increase in the field strength. The same behaviour takes 
also place for K > 1 when the perturbative analysis may be not trusty. But, as 
we have determined, the picture of the symmetry behaviour is only quantitatively 
changed for heavy scalar particles: the first order phase transition in weak fields 
becomes the second order one for strong fields. These circumstances convince us 
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that the assumption of Ref. that the hypermagnetic field makes the weak first- 
order phase transition strong enough is not proved by the detailed calculations. 

In Table 2 we adduce the results for the phase transition obtained within the 
asymptotic EP. 



h 


K 


T c (GeV) 




<pc{h,T c ) 


R 


M w {h,T c ) 


0.01 


0.85 


104.678 


0.097 


0.3114 


0.7319 


0.3358 


0.01 


1.25 


119.569 


0.098 


0.3146 


0.6473 


0.3441 


0.01 


2 


142.01 


0.102 


0.3194 


0.5532 


0.3563 


0.05 


0.85 


105.449 


0.085 


0.2915 


0.6801 


0.3050 


0.05 


1.25 


120.241 


0.094 


0.3066 


0.6765 


0.3319 


0.05 


2 


142.559 


0.102 


0.3194 


0.5511 


0.3611 


0.1 


0.85 


second 


order 


phase 


transition 




0.1 


1.25 


121.616 


0.062 


0.2490 


0.5037 


0.2809 


0.1 


2 


143.539 


0.088 


0.2966 


0.5084 


0.3420 



Table 2. 

These data, as previous ones, show that the ratio R is less then unit and the 
baryogenesis condition does not hold. 

Now, let us compare the results in Tables one and two. As one can see, even 
for h = 0.01 they differ considerably for all the parameters exept the critical 
temperatures. The first order phase transition for K = 0.85 determined by both 
of the potentials possesses the same characteristics. But this is not the case for 
K = 1.25, 2 when the phase transition described by the exact EP is weaker first 
order. For the field strength h = 0.1 the exact EP predicts the weak first-order 
phase transition for all K whereas the asymptotic one fixes the second order 
phase transition for K = 0.85. For K = 1.25, 2 the jump of the order parameter 
determined by the asymptotic EP is twice larger then for the exact EP case. If 
the fields are stronger then 0.2 - 05 -10 24 G both EP predict the second order 
phase transition. Thus, we conclude that asymptotic EP predicts conversion of 
the first order phase transition to the second order one for weaker fields whereas 
the first order phase trasnsition described by it is stronger as copmare to the 
transition derived within the exact EP. 

9. SYMMETRY BEHAVIOUR IN STRONG MAGNETIC 
FIELD 

As we have described qualitatively in sect. 2, the case of external magnetic 
field somehow differs from the hypermagnetic one. Nevertheless, a nubmer of 
previous results dealing mainly with the form of the EP curve in the broken 
phase are relevant, since in this phase the unscreened constituent of the external 
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hypermagnetic field coinsides with the magnetic field. The condition fixing the 
transition temperature is 

Vrestored{H,T c ,0) = Vb ro ken(H , T c , (j) c (H , T c )) . (71) 

The transition happens when the depth of the minima located at the begining, 
4> c = 0, and at (fi c ^ is the same. 

Below, as before, we present the results of numeric investigations of the phase 
transition obtained for the exact EP and for the high temperature limit of it. In 
Table 3 we show the characteristics of the first order phase transition determined 
within the exact EP. 



h 


K 


T c (GeV) 




4> c (h,T c ) 


R 


M w (h,T c ) 


0.01 


0.85 


105.18 


0.107 


0.3271 


0.76506 


0.3504 


0.01 


1.25 


120.77 


0.045 


0.2121 


0.4321 


0.2541 


0.01 


2 


143.96 


0.016 


0.1265 


0.2162 


0.2031 


0.1 


0.85 


106.35 


0.098 


0.3130 


0.7241 


0.2835 


0.1 


1.25 


121.92 


0.021 


0.1449 


0.292 


0.1243 


0.1 


2 


146.14 


0.003 


0.0545 


0.0917 


0.1272 


0.5 


0.85 


108.23 


0.092 


0.3033 


0.6894 


0.4696 i 
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1.25 


second 
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phase 


transition 
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second 
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phase 


transition 





Table 3. 

The second column shows the values of the parameter K = m 2 H /M^ } corre- 
sponding to the Higgs boson masses 75 GeV, 90 GeV and 115 GeV, respectively. 
All the parameters are as in section 8. 

It is seen, an increase in h makes the first order phase transition weaker, as 
in the case of hypermagnetic field. For strong fields H > 10 22 — 10 23 G the 



baryogenesis condition R > 1.2 — 1.5 [13] is not satisfied. 



Now, let us show the results for the asymptotic EP. They are gethered in 
Table 4. 
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h 


K 


T c (GeV) 


4> c (h,T c ) 


R 


M w {h,T c ) 


0.01 


0.85 


106.321 


0.0316 


0.0732 


0.1310 


0.01 


1.25 


122.022 


0.0548 


0.1104 


0.1517 


0.01 


2 


145.959 


0.0616 


0.11 


0.1640 


0.1 


0.85 


107.267 


0.0316 


0.073 


0.1319 


0.1 


1.25 


second order 


phase transition 






0.1 


2 


second order 


phase transition 







Table 4. 

First of all notice, the baryogenesis condition R > 1.2 is not satisfied for 
neither the exact EP nor the asymptotic one. The R values are small for for all 
K considered. Since the external field is accounted of exactly through the Green 
functions the effects of the field influence are also exactly reproduced. Thus, 
our analysis has shown that in the case of external magnetic field baryogenesis 
can not be generated in the SM. The comparison of the results for the exact 
EP (Table 3) and the asymptotic EP (Table 4) shows that in the latter case the 
second order phase transition is predicted for more weak external fields. 

The main results of this and previous section can be summarized as follows: 
the external either the magnetic or the hypermagnetic field can not produce the 
strong first order EW phase transition. Baryogenesis does not survive in the SM 
due to these external conditions. 

In the next two sections we shall investigate the problems dealing with the 
spontaneous magnetization of the vacuum at high temperature and stabilization 
of this state due to radiation corrections. We shall show that strong magnetic 
fields H ~ g^ 3 T 2 have to be spontaneously created. 

10. RADIATION SPECTRA OF W-BOSON IN STRONG 
MAGNETIC FIELD AT HIGH TEMPERATURE 

VT-boson mass operator in a constant magnetic field H at zero temperature has 
been calculated and investigated in Refs. |)9| , [|7(J . In particular, it gave possibil- 



ity to clarify the role of the radiation corrections in the problem of stabilization of 



the jy-boson spectrum (see survey |70| ). The temperature dependent radiation 



corrections to the M^-boson spectrum has been studied in Ref. |}6|. The longi 



tudinal components of gauge fields acquire the temperature masses ~ gT [|15 
The tachyonic mode is the transversal state excited by a spin interaction. Its 
tepmperature mass can be calculated as the average value of the mass operator 
in the ground state of the spectrum. 

To incorporate temperature the imaginary time formalism will be used. As 



in the case of zero temperature [B9|, the Schwinger operator method and s- 



representation will be applied. In general, after summation over discrete imagi- 



32 



nary frequencies this method becomes not applicable. However, it remains prac- 
tically unchanged in the limit of high temperature when only the static modes 
I = contribute [[HJ. This approximation is sufficient to investigate the role of 
the daisy diagrams and will be used in what follows. 

We calculate the average value of the W-bosoia mass operator in the states 
of jy-boson spectrum in a magnetic field | n, a >, where n, a are the Landau 
level number and the spin projection variable, respectively. These functions, < 
M(H,T,n,a) >, give the radiation temperature dependent masses of the states. 
The effective mass squared is M 2 (H, T) = M 2 — eH + Re < M(H, T, n = 0, a = 
1) >. If this value is positive, the spectrum is stabilized by radiation corrections. 

In the present section we consider a simplified model of electroweak interac- 
tions (the boson part of it) based on the spontaneous breaking of SU(2) — > U(l) 
gauge group. The Lagrangian is 

L = -\F^ + \(D^r + ^-^\ (72) 

where x 2 = x a x a ,a = 1,2,3, F£ v = d^A a u - d u A^ + ge abc A b ^A c u , D ab = d^5 ab + 
ge abc A c . This model is described, for example, in Ref. ]7D[. After the sponta- 
neous symmetry breaking the charged components = ^(A^ ± iA 2 ^) asquir 

the masses M = g<f) c , where <\> c =\ m \ /\f\ is the vacuum value of the scalar 
field, and the component A^ = A^ remains massless and is identified with the 
electromagnetic potential. We also identify the gauge coupling constant g with 
electric charge: g = e. In fact, this model is the mass regularization of the Yang- 
Mills theory. In the limit M — > all the corresponding results are reproduced. 
The results for Salam- Weinberg theory can also be obtained by the simple sub- 
stitutions, 7-contributions — > Z-contributions with the corresponding vertex and 
propagator factors (see for details Ref. |ffO| ). 

To investigate the problem we are interested in, let us, as before, direct the 
external magnetic field along the third axis H = if 3. The corresponding potential 
is chosen in the form A e * 1 = (0, 0, Hx 1 , 0), H = const. To quantize the fields the 
following gauge fixing condition is used: 

d^W^ - ieAfW^ - M<j)± = 0, (73) 

where ^ = ■^(4 )1 =t ^0 2 ) are the charged Goldstone fields. 

The V^-boson mass operator in one-loop order is given by a standard series 



of diagrams ||70|| , |69| . It can be written in the form 



&4 



where (3 = 1/T, k 4 = 2?rZ//3, I = 0, ±1, ±2, operator P„ = id^ + eA^K The 
first term in the Eq. ( JT^ j describes the contribution of virtual neutral scalar 
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particles and the second one gives that of virtual gauge fields (photons, W- and 
Goldstone bosons). 

We restrict ourselves by investigation of the high temperature limit, that cor- 
responds to the static mode 1 = 0. In this case the standard calculation procedure 
developed at zero temperature can be straightforwardly applied. Details of these 



calculations can be found in Refs. |j9) . J7(J , |66| . Here, we note the points specific 
for /c 4 = case. For brevity, in what follows we will consider the part M^(k, P). 

To integrate over d 3 k in Eq. ([F4]) one has to introduce the s-representation 
for each propagator entering and present the product of propagators as follows 

1 oo 

D^(k)G w (P -k) = - J duf dsse isH , (75) 



where " Hamiltonian" is, 

H = (1 - u)k 2 + u(P - k) 2 - u(M 2 + 2ieF) - (1 - u)m 2 , (76) 

and m is the mass of the neutral scalar field. The ^-dependent part of the Hamil- 
tonian is zero. Then, the integration is carried out by means of the procedure 
developed in Ref. pTf. The difference between the finite and zero temperature 



cases consists in the dimension of the corresponding integrals. This is reflected in 
the power of s appeared after integration. Three dimensional integral (at T ^ 
) gives the factor s™ 3 / 2 whereas at T = one obtains s~ 2 ||61||. 



To find the energy of the states owing to the radiation corrections one has to 
define the jy-boson mass-shell in the field at high temperature. In the static case 
it is described by the equations 



[ ( P 2 + M 2 )5 i3 + 2ieFZ* ] W~ = 0, (77) 



pext 

ij 1 ./ 

pwr = o, <tr = 0, 



where j = 1,2,3 and the product P4H4 = bacause for static modes p 4 = 0. 
The states are normalized by the condition 

< n, a I n', a' >= 5 n ^ 8 a y, (78) 

n, n' = 0, 1, ... and a, a' = 0, ±1. 

The average value of the mass operator in these states can be written in the 
form 

1 00 

< u > = [duf —j= [ettiif A] -y 2 e ~ xuM2 / eH (79) 

v r v 

exp [-(2n + l)(p - x{\ - u) - 2u(l - u)]M(x,u), 
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where M(x, u) =< n, a \ Mij | n, a > and 

(1 — u)shx , . 

tanhp = - \ , ; — — 80 

(1 — u)chx + usnx/x 

A = (i_ M )2 + 2n(l-n)^ + n 2 (— ) 2 . 

2x x 

This very complicate expression can be investigated for different limits of interest. 
In particular, for high temperatures and strong fields, |j| >> 1, << 1 we 
obtain (in the reference frame = 0): 

e 2 

< n,<7 = +1 | M 9auge | n,a = +1 > = — (e#) 1/2 T[12,33 + 4n (81) 

An 

+ i(3 + 6n)], 
e 2 

< n,a = +1 | M 4, | n,(7 = +1 > = — (ei/) 1/2 r[14, 63 + An + z(7 + 6n)], 

i?e<n,(T = -l|M sm9e |n, ( T-l> = — (e#) 1/2 T(ll, 44 + 4n), 

An 

e 2 

i?e < n,cr = I M 9auge I n,a = > = — (e#) 1/2 T(15, 44 + 4n). 

An 

Above formulae give a general picture on the behaviour of the radiation energies 
for different spin states. As it is seen, the real part of < M > is positive in the 
ground and excited states. It acts to stabilize the tree spectrum. The imaginary 
part describes the decay of the states due to transitions to the ones having lower 
energies. 

Now, let us consider the effective W^-boson mass squared: 

M 2 (H,T) = M 2 — eH (82) 
+ Re < n = 0, a = +1 | M gauge + M 4, | n = 0, a = +1 > 

e 2 

= M 2 -eH + 26,96— (eH) l/2 T. 

An 

This value is positive for sufficiently high temperatures. Hence, one can conclude 
that radiation corrections in the field stabilize the vacuum at high temperatures. 
The temperature mass of the transversal modes dependens on the field and equals 
to zero when H = 0, as it should be \TE\. 

Obtained results are of interest for cosmology. Namely, if at the EW phase 
transition epoch the magnetic field was present, the radiation mass of W^-bosons 
could serve as the dynamic mechanism of the vacuum stabilization in both the 
broken and the restored phases. It was discussed in detail in the previous sections. 
To derive a consisten picture one has to consider a vacuum magnetization at high 
temperature with the correlation corrections taken into consideration. 
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11. THE SPONTANEOUS VACUUM MAGNETIZA- 
TION AT HIGH TEMPERATURE 

The generation of magnetic fields in nonabelian gauge theories at finite temper- 
ature is of great importance for particle physiscs and cosmology. Its positive 
solution, in particular, will give a theoretical basis for investigations of the QCD 
vacuum at high temperature and the primordial magnetic fields in the early uni- 
verse ||20|| , ||36|| . In literature different mechanisms of producing the fields are 
discussed (see recent papers [|J, |L3| and references theirin). In this section we 



investigate in more detail one of possibilities - the spontaneous magnetization 
of the vacuum of nonabelian gauge fields at finite temperature. This problem 
has been studied in one-loop order in Refs. [20], [17|, where the creation of 



the vacuum field was derived. In Ref. |L7| a number of correlation corrections 
has been taken into account, but the polarization functions of gauge fields were 
not calculated and therefore a trusty conclusion about the phenomenon was not 
obtained. In Ref. [|ll[] that has been considered with all necessary daisy diagrams 



included. Below, we shall follow this paper. 

For simplicity, we investigate the vacuum magnetization at finite temperature 
within SU(2) gluodynamics. Considering the abelian covariatly constant chro- 
momagnetic field H a = S a3 H = const and finite temperature as a background, 
we calculate the EP containing the one-loop and daisy diagram contributions of 
the neutral and charged gluon fields. The field dependent Debye mass of neutral 
gluons can be computed from a special type EP. To find the one of the charged 
gluons the high temperature limits of the gluon polarization functions at the 
background are also calculated. It will be shown that in the adopted approxima- 
tion the Savvidy level with the field strenght (gH)^ 2 ~ <? 4//3 T is generated. This 
is strong field. Although it is screened at distances I > (g 2 T) _1 by the gluon 
magnetic mass, the spectrum of charged particles, being formed at Larmor's ra- 
dius r ~ (gH)c 1 ^ 2 ~ ((7 4 / 3 T) _1 , is located inside this domain, r << l m , for small 
9- 

Most of results for the Yang-Mills (YM) theory can be obtained without actual 
calculations by making use of the results of the previous chapters. Really, the 
model considered in sect. 9 is a mass regularization of the YM theory. The 
external magnetic field just coinsides with the abelian chromomagnetic field of 
interest. Therefore, the results for the latter one can be obtained by setting 
the mass of the charged gauge fields to zero and subtracting the contributions 
of the longitudinal spin projection of the massive charged gauge fields and the 
Goldstone fields. These transformations can be made freely. For instance, to get 
the contribution of charged gluons to the one-loop EP it is sufficient to substitute 
the factor 3 by the factor 2 in the part of Eq. (|20|) describing the contribution of 
massless (M = 0) charged spin zero particles. The same has to be done to find 
either the Debye mass of neutral gluons (in Eq. (|48|)) or the temperature masses 
of charged gluons (sect. 9), etc. Below, for convenience of further account, we 
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adduce the necessary expressions [TTI . 



The high temperature limit of the one-loop EP is 

- i^^ + 0{g 2 H 2 ). (83) 

Notice the cancellation of if -dependent logarithms entering the vacuum and the 
statistical parts, /i is a normalization point. 

11.1 Daisy diagrams of the unstable mode and the neutral 
gluons 

The contribution of daisy diagrams with the unstable mode is given by the ex- 



pression |l I 



VunstaUe = ^[II(ff , T, U = 0, a = +1) - gH] 1 ' 2 + l iiE^L . (84 ) 
Z7T Z7T 



From Eqs. (^J) and ( jS4|) it is seen that the imaginary terms are cancelled out in 
the total. The final EP is real if the condition H unstc ^i e {H ,T) > gH holds. 

This expression must be supplemented by the term describing the contribution 
of the neutral gluon fields. In one-loop order it gives a trivial if -independent 
constant which can be omitted. However, these fields are long-range states and 
they do give if -dependent EP through the correlation corrections including the 
mass term II (if, T). Corresponding part of the EP is described by the expression 



which can be recognized from the case H = |46[|, [[To" 



V rinq = —U°(H,T)T 2 l —[Y[\H,T)f 2 

9 24 v ; 12tt/3 l v n 

+ 32,2 [Mp^ + j-q, (85) 

where n°(if, T) = IIq (A: = 0,H,T) is the zero-zero component of the neutral 
gluon polarization operator calculated in the external field at finite temperature 
and taken at zero momentum, C is Euler's constant. The first term in Eq. ( j35|) 
has order ~ g 2 in coupling constant, the second term is of order ~ g 3 and the 
last one - ~ g A . Restricting ourselves by order ~ g 3 , it will be omitted in what 
follows. As usually, for II (if, T) the high temperature limit of the function has to 
be used. The expression (^) needs in one additional comment. If one compares 
it with Eq. (|35D , the difference will be in the extra term ~ II(T, if )T 2 . We have 
maintained this next-to-leading term because it is of great importance for prob- 
lems under consideration. In studying of symmetry behaviour only the leading 
(^-independent) terms of the polarization functions were taking into account. 
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The mass squared m 2 D neutra i = Uoo(H,T,p = 0), that is the Debye mass of 
neutral gluons, reads 



m 2 D = \g 2 T 2 - - -LfofO + 0((gH) 2 /T 2 ). (86) 

3 7T AlT Z 

Here, the first term is the well known temperature mass squared and other ones 
give the field-dependent contributions. They have negative signs that is impor- 
tant for what follows. Sabstituting expression (|86"D into equation (|S5|) we obtain 
the correlation corrections due to neutral gluons. 

The calculations, as we have described in sect 10., result in the high temper- 
ature limits of U unstab i e (H,T) 

UunstabUH, T) = <n = 0,a = +l\ W^ d \ n = 0, a = +1 > 

a 2 

= 12.33^(gH) 1/2 T, (87) 
and of excited states U(n ^ 0, cr), 

a 2 

ReU(p 4 = 0,n,p 3 = 0,H,T,a = +1) = ^-{gH) 1/2 T {12. 33 + An), (88) 

Att 

a 2 

ReU(p 4 = 0,n,p 3 = 0,H,T,a = -1) = gH) 1/2 T( 11. AA + An), 

Att 

where the average values of the polarization operator in the corresponding states 
of the spectrum (|9l|) are computed. These formulae have been obtained in the 
high temperature limit gH/T 2 << 1. The operator contains also an imaginary 
part which describes the decay of the states. But for the problem under con- 
sideration only the real part is needed because it is responsible for the radiation 
masses of particles. 

Let us note the most important features of the expressions (|8~7D , ([38l). It is 
seen, at H = no screening magnetic mass is produced in one-loop order. Second, 
the mass squared of the modes are positive and act to stabilize the spectrum of 
charged gluons at high temperatures. Therefore, in the nonzero chromomagnetic 
field at finite temperature the charged transversal gluons become massive. 



11.2 Contribution of the daisy diagrams of charged gluons 

To obtain the correlation corrections due to the logitudinal charged gluons the 
zero-zero component of the polarization operator has to be calculated. The Debye 
mass of charged gluons is found to be 

n 00 (A;4 = 0, h = 0, H, T) = 2 -g 2 T 2 + ^-{gHf^T^ + An), (89) 

where again only the real part is adduced. Hence, the next-to-leading terms are 
the growing positive functions of n. 
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Now, let us turn to the generalized EP 
^ = f E T| E MW + <L P3 + n(n,T, H,a))] (90) 

Z7r i = _ 00 _^ o Z7r n=0 ,cr=±l 

written as the sum of energies of the charged gluon field modes in the external 
chromomagnetic field ( coi = - discrete imaginary energies) 

el a =pj + (2n + l- 2a)gH + U(H, n, a, T) (91) 

and including the temperature masses U(H,n,a,T) of the modes which are also 
dependent on H, the level number, n, and the spin projection a. In general, the 
incorporation of the polarization functions into EP may result in wrong combi- 
natoric factors for the two-loop diagrams. In our case, to be sure in the obtaining 
results the following procedure is applied. We subtract the term 

V. = Y L Y,fp- E D (p 3 ,H,T)U(H,n,T,a) (92) 

l=— oo_ OQ n=0,cr=±l 

from Eq. (|90|) , that separates the contributions of the two-loop diagrams of 
charged gluons. One freely can check that the two-loop diagrams containing the 
one of loops with the neutral gluon lines can be accounted of within the general- 
ized EP including the insertion of the neutral gluon polarization operator. After 
that the only two-loop diagram (two touch circles) allowing for the self-interaction 
of charged gluon fields in the vacuum remains to be computed separately. 

Then, by substituting the expressions (|57D , fl88D , (|8^) into Eq. ( p0|) and 
integrating over momentum and calculating the sums in n, we obtain the daisy 
diagram contribution of charged gluon fields. The result can be expressed in 
terms of the generalized (^-function and looks as follows, 

Kt 9 = 9 ^{\f^D\a-\^ + )+a-\^-)+K{-\ a D )] 

+ (U(H,T,n = 0,cr = +1) 1/2 }, (93) 

where the first term in the first squared brackets corresponds to the spin projec- 
tion a = +1, the second term - a = — 1 and the last one describes the part due 
to longitudinal charged gluons. The terms in the second squared brackets give 
the independent of H(H, T) part of eq. ([K]). The last term in the curly brackets 
is due to radiation mass of the unstable mode. The notations are introduced: 

yj^n U rL ^2iT\y rL > 1 ' "~ 2 ^ 4tt 2gH D ' "+ 2 ^ 4tt 2gH D d,UU 

a D = (\g 2 T 2 + ^{gHf/ 2 T + gH)/2gH D . This expression is real for sufficiently 
high temperatures. 
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Substituting the expression flSBD in eq. fl55| ) and gethering all other contribu- 
tions, we obtain the consistent EP. 

To calculate the two-loop vacuum diagram describing the self-interaction of 
charged gluons the standard procedure can be applied. The computation details 



are given in Ref. ]TTJ. Reffering the readers to this paper, we write down here 



VS\H,T) IT _ = -^-{gHf/V + 0{gHT 2 ). (94) 



the resulting expression for the high temperature limit 

f 
Air 

The high temperature limit of the term subtracted from the generalized EP fl92]) 
is 

V s (H,T) lT _ = ^-C(^,-^)(^) 1/2 T 3 + 0(^T 2 ). (95) 

It also gives the negative contribution to the leading terms of the asymptotic 
expansion since — = 0.8093. It worth to mention that these are the lon- 

(2) 

gitudinal modes that determine the high temperature behaviour of V^ h (H,T). 
Having obtained the two-loop corrections to the EP, one can investigate the spon- 
taneous vacuum magnetization at high temperature. 



11.3 Vacuum magnetization and the stability problem 

The derived EP is expressed through the well studied special functions. Therefore, 
it can easily be investigated numerically for any range of parameters entering. 
As usually, it is convenient to introduce the dimensionless variables: the field 
4> = (gHy/ 2 /T and the EP v((f),g) = V(H,T) /T 4 . The vacuum magnetization 
at high temperatures, T >> (gH) 1 ^ 2 ,^ — > 0, will be investigated within the 
followin limiting form of the EP, 



2g 2 48 7T 2 "> 2/ 3 7T 48tt~ 



3 3 ; vr 6+f 4vr Qy/2- 
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0,80930, (96) 



where other ~ g 3 terms are omitted. The logarithmic term is signalling the 
asymptotic freedom of g 2 {T) at high temperatures It includes explicitely 



the dependence on the scale parameter /x. Other terms present, respectively, the 
high temperature limits of the one-loop EP, the neutral gluon and the charged 
gluon daisyes and the two-loop diagram of charged gluons. To obtain the term 



ch 

1 



due to V™ 9 the asymptotic expression for Zeta-function (72 



C(-|, Mla^oc = -| fl D 2 + ~ 4^ a D /2 + °( a D /2 ) ( 97 ) 
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was used. Zeta-functions with a+,a_ do not contribute in leading order. Since 
we are searching for the fields of order larger then g 2 , we can omit the term g 3 
and obtain for the condensed field 

(9H)l /2 = % 9 4/3 T. (98) 

Thus, we come to the result that the ferromagnetic vacuum state indeed exists 
at high temperatures. The correlation corrections increase the field strength as 
compare to the one- loop value ( (gH) l J 2 ~ g 2 T). 

Let us discuss the stability of the condensed field. The second derivative 
of the EP is positive for H c that means we have the minimum. The field is 
not changing in the direction a = 3 of the isotopic space. To check that this 
is indeed the case for the perpendicular directions = 1,2 or a = responsi- 
ble for excitation of charged fields W ± , one has to calculate the effective mass 
squared M 2 (H C ,T). First we consider the one-loop case. Sabstituting the value 
(gH^) 1 / 2 = (g 2 /2n)T in the one-loop polarization function, we find that the 
effective mass squared, M 2 (H^\T) = U^(H^\n = 0, o = +1) - gH^ > 0, 
is positive. Thus, the vacuum stabilization is observed in this consisten calcula- 
tion. However, if one checks whether the one-loop gluon radiation mass stabilize 
the true vacuum magnetic field and substitutes the value (gif c ) 1//2 ~ g 4 ^ 3 T, the 
negative value of M 2 (H C ,T) will be obtained. The one-loop mass does not sta- 
bilize the spectrum and, hence, vacuum. Nevertheless, the one-loop result makes 
hopeful the idea to have the stable vacuum due to radiation corrections to the 
charged gluon spectrum. Naturally, to investigate this possibility the gluon po- 
larization operator with the correlation corrections included should be calculated. 
This problem requirs an an additional investigation. Other interesting possibility 
is the formation at high temperatures of the gluon electrostatic potential, so- 
called A Q condensate (see survey fl68"l), which also acts as a stabilizing factor 



To realize the latter scenario consistently the simultaneous spontaneous gener- 
ation of both the Aq condensate and the magnetic field should be investigated. 
If again the homogeneous vacuum field will be found to be unstable with these 
improvements made, the inhomogeneous fields of the lattice type discussed in 



Refs. |p6| , |)43l , [pq| can be created. Really, since the condensed magnetic field 
is strong at high temperatures, the lattice structures having the cells of order 
~ l/(g 4 / 3 T) << \/{g 2 T) are located inside the domain where the fields are not 
screened by the gluon magnetic mass. 

The above calculations have unambiguously determined the possibility of the 
vacuum magnetization at high temperature, although a number of questions con- 
cerning the vacuum stabilization and structure has to be investigated in order to 
derive a final picture. This could result in the presence of strong magnetic fields 
in the hot universe. 
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12. COMPARISON WITH OTHER APPROACHES 



In this section we are going to compare our results with that of other investiga- 
tions. That refer to either the phase transition or the generation of the magnetic 
field at high temperature. We begin with discussion of the EW phase transition 
in strong hypermagnetic fields. 

First let us discuss the results of Refs. 0, || where it was concluded that 
strong external hypermagnetic field generates the sufficiently strong first order 
phase transition and baryogenesis survives in the SM. In these papers the inflence 
of the field on the vacuum has been allowed for in tree approximation that gives 
the qualitative estimate of the effect. In Ref. || the temperature dependent part 
of the EP has been included in one-loop order whereas the field dependence was 
skipped. Naturally, further studying of the phenomenon should include the field 
dependent radiation and correlation corrections due to fermions and bosons. This 
is the problem that we have addressed to in the present investigation. The main 
idea was to determine the form of the EP curve in the broken phase and find 
the range of the parameters Hy, K when the EW phase transition of first order 
is happened. To elaborate that the consistent EP including the one-loop and 
daisy diagrams of all the fundamental particles has been calculated. As we have 
discovered, the role of fermions is crucial in a vacuum dynamics in strong fields 
at high temperature. They essentially affect the structure of the broken phase 
making the EW phase transition weaker as compare to the tree level results. 
The external field has been accounted of exactly through Green's functions. The 
minimum of the EP was found to be real at sufficently high temperatures when 
the first order phase transition happens. This important property was established 
when the daisy diagrams of the tachyonic mode have been included. As a result, 
no conditions for W- and Z-boson condensates are observed at high temperatures 
and the external field strengths corresponding to the first order phase transition. 
The condensates could be generated for stronger fields when the phase transition 
becomes of second order. But, this is not of interest here since we are looking for 
conditions when baryogenesis can survive. 

In Refs. ||, |)7]] the EW phase transition in the hypermagnetic field has been 



investigated by means of a general method developed in Ref. [FFl] (see also survey 



H|) which combines the perturbative computations and the lattice simulations. 
The results obtained therein have supported the main conclusions of Refs. H, 
discussed above. In more detail, it was found for the fields Hy < 0.3T 2 the 
first order phase transition becomes stronger, but it still turns into a crossover 
for masses ran > 80 GeV. For stronger fields, a mixed phase analogous to a first 
type superconductor with a single macroscopic tube of symmetric phase, parallel 
to Hy, penetrating through the broken phase, has been observed. 

As this aproach is concerned, we note that because of peculiarities of the 
calculation procedure adopted in Ref. || all fermions except t-quarks are decou- 
pled as nonstatic modes and therefore the field dependent fermion contributions 
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as well as the correlation corrections could not be accounted of. That is why 
these calculations also do not reproduce correctly symmetry behaviour for strong 
external fields. 

In the present approach, the value of the mass uiu = 75 GeV ( = 0.85) is close 
to ran = 80 GeV discussed in Ref. ||. This is important for us because the first 
order phase transition is controlled by perturbative method used. But we could 
not observe the mixed phase. We have seen the crossover (or second order phase 
transition) for Hy = 10 23 G and K = 0.85 determined within the asymptotic EP. 
The exact EP in this case predicts the weak first-order phase transition. Since 
the EP is real in the minimum, no conditions for the vortex-like phase exist. 

The Higgs boson mass values considered correspond to the cases when per- 
turbative results are reliable (K = 0.85) and may be not trusty (K = 1.25 , 
2). However, since the external field is taken into account exactly its effects 
are correctly reproduced. As we have seen, an increase in Hy makes the EW 
phase transition of second order for field strengths Hy ~ 0.5 • 10 24 G for all the 
values of K investigated. For weaker fields the phase transition is of first order 
but the ratio R = <f) c (H,T c )/T c is less then unit, that is unsufficient to generate 
baryogenesis. 

Let us remind the situation with the magnetic field when the magnetic mass of 
order m mag ~ g 2 T is taken into consideration. This mass screens the nonabelian 
component H5 a3 at distancies I > l m ~ (g 2 T)~ l but inside the space domain 
I < l m it may exist and affect all the processes at high temperatures. The latter 
fact has not been realized in a number of investigations [33], [|35| where (as in 



Ref. |2(J ) the field strength generated at finite temperature has been erroneously 
estimated as coinsiding with that at zero temperature. Our estimate of the field 
strength at high temparature makes the investigation of the EW phase transition 
in strong magnetic fields resonable. But, as we have seen in sect. 9, in this case 
it also is impossible to generate the strong first order phase transition. 

Now, let us consider in more detail the problems of the magnetic field gen- 
eration and stabilization at high temperature. In the present survey we have 
investigated in detail the Savvidy mechanism. The ways of the vacuum stabi- 
lization have also been discussed. In this scenario, usual magnetic field should 
be treated as the projection of the chromomagnetic field created in a nonabelian 
gauge theory. Results of our investigation disagree with that of Refs. [|34| , |35 



where it is claimed that spontaneous magnetization does not hold at high temper- 
ature. To clarify the origin of the discrepancy let us repeat the main statements 
of Refs. p4] |, p5| : 1) The field strength generated at finite temperature g H ~ A 2 



where A 2 = /i 2 exp(— u^z^ ) , coinsides with that at zero temperature. This is 
much less than the magnetic mass squared ~ (g 2 T) 2 . 2) Since the spectrum of 
charged particles in the magnetic field is formed at the Larmor radius scale, the 
weak long range fields are not produced being screened by the mass at distances 
I — ^-/g 2 T- The error is the assumption that the field strength generated in the 
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vaccum is not changed when the temperature is switched on. As we have seen 
in sect. 11, strong abelian colour magnetic fields of order (gH) 1 ^ 2 ~ g 2 T (in 
one-loop approximation |T7[ , Jl9[| ) or (gH) 1 ^ 2 ~ g 4 ^T (Tl| (when higher order 
corrections are included) is spontaneously generated at hight temperture. 

We also have observed in the consistent calculation that the Sawidy state 
is stable at high temperature when the one-loop EP and the one-loop gluon 
polarization operator are taken into account. 



13. DISCUSSION 

As it was realized recently, in the SM the usual scenario of baryogenesis can not be 



established without external fields By analogy to superconductivity it was 
assumed that strong hypermagnetic field is able to generate the strong first order 
phase transition for the values of mass m# permitted by experiment ||. Different 
mechanisms to produce magnetic fields have been proposed in literature. They 
can be devided in two groups: 1) generation of fields at EW phase transition p2|], 
fg, [g|, f2§; 2) generation of fields beyond the SM scale (for instance, GUT 
theories) @], A special interest was in strong hypermagnetic field which due 
to its abalian nature is not screened at high temperature whereas the nonabelian 
component of the usual magnetic field is screened at distances / > l/g 2 T by the 
gauge field magnetic mass. The influence of the magnetic fields generated at the 
EW scale on the first order phase transition was investigated in Ref. |27|. It 



has been shown that if the field strength is stronger then Hq = 10 24 G the phase 
transition is delayed. Of cause, it is relevant if the first order phase transition 
is realized. But this is not the case for the SM, where without external fields 
for the mass itlh > 60 — 70 GeV the second order phase transition is predicted. 
Naturally, internal forses are not able to change the kind of the phase transition. 

The investigations carried out in the present paper are refered to the second 
possibility assuming that the fields have been generated before the EW phase 
transition at a GUT scale. The dynamics of the hypermagnetic field generation 
has not been discussed here, although the corresponding machanism was proposed 
in Ref. 0]. We just assumed that the field present in the early universe. As the 
magnetic field is concerned, in sects. 10, 11 we have investigated in detail its 
generation and stabilization at high temperature. 

The results on the hypermagnetic field influence can be summarised as follows. 
In contrast to the conclusions of Refs. |J, || claiming that strong hypermagnetic 
field generates the strong first order EW phase transition, we observed an opposite 
effect. As we have seen in our numerical computations, the weak first-order EW 
phase thansition bacomes of second order in strong fields. The origin of the 
discrepancy is the following. These authors have not determined correctly the 
form of the EP curve at the transition temperature since the contributions of 
fermions and the correlation corrections in the fields were skipped. In fact, our 
analysis has completed the investigations not only for the hypermagnetic but 
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also for the magnetic field. We have discovered that the role of fermion radiation 
corrections is of great importance in the phenomenon investigated. To better 
understand the role of fermions in symmetry behaviour let us adduce two terms 
of the asymptotic expansion of the EP in the limit of T — > oo, H —>■ oo. The first 
one is the term ~ H 2 log^j. Due to this term the light fermions are dominating 
at high temperature. The second term can be derived from the expansion of the 
zero temperature part Eq. fl32|). This expression side by side with the leading 
term ~ H 2 log^, which due to a "dimension parameter trading" is replaced by 

the above written term, contains the subleading one ~ —eHm'jlog^j (for details 



see Ref. [2I|). This term acts to make "heavier" the Higgs particles in the field. 
As a result, the second order temperature phase transition is stimulated due to 
strong fields. We also have investigated the influence of different parts of the EP 
on symmetry behaviour. It was discovered that the change of the kind of the 
EW phase transition with increase in H is due to the fermion temperature part 
of the EP. 

The stability problem of the vacuum at high temperature and strong fields 
has also been investigated. Our EP (its minima) is real at the first order phase 
transition that is important for the reliability of the results obtained. The lat- 
ter property is insured by the imaginary terms of daisy diagrams cancelling the 
corresponding ones of the one-loop EP. Thus, the total EP is sutable to correctly 
describe the phase transition. As we have found, the effective W-boson mass 
M W (H,T) is real in the local minimum, therefore no conditions for the forma- 
tion of the vortex-like structure due to W- and Z-boson condensates exist in the 
broken phase. Probably, this is the reason why the condensates have not been 
observed in Ref. ||, although it is difficult to check that strightforwardly. 

Furthermore, we have also shown that spontaneous vacuum magnetization in- 
deed happens at high temperature giving the reliable mechanism of the magnetic 
field generation at a GUT scale. In principle, the field could be spontaneously cre- 
ated at the EW scale. However, in this case it is incorrect to treat the EW phase 
transition as the external field problem. It is necessary to apply consideration as 



in Ref. 27 



On the base of our analysis one has to conclude that baryogenesis can not 
survive in the SM if smooth external magnetic fields generated beyond the EW 
scale are included as environment. 

As a final remark, we would like to stress once again that light fermions act to 
turn the first order phase transition into the second order one with increase in the 
field strength, independently of the values of the Higgs boson mass considered. 
This fact may be one of the reasons why it is of interest to carry out a similar 
investigation for the minimal supersymmertic SM or other super symmetric ex- 
tensions of it, because the fermion and the boson sectors enter such models on 
an equal footing. 
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APPENDIX 

General method of calculation of high temperature asymptotics has been devel- 
oped in Ref. fl64"| . To outline the procedure used in the present paper let us 
consider the fermion contribution as an example: 



n=l 

At first, let us make Mellin's transformation of K\ r \{nuj) with respect to parameter 



n 
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OO I I I I 

K M (a) = J K^f-'dt = ^ s 2 s " 2 r(| - y) r (i + y)' 



where T(x) is T-function, and substitute it in Eq. (A.l). Then we find 

/^EM)" / dsl! --(>r(£-M) r (£ + M), (A2) 

Tl — 1 s~i 

— ioo 

where C r > \r\ and the integral is calculated along the straight line parallel to 
the imaginary axis. If one shifts the integration variable s — > s — r = s', the 
integration contour moves to the right, C<i r — ► C 2r , for negative r, and to the left, 
C r — > C' , in the case of positive r: 2r < C' 2r <2r — 1; < Cq < 1. Changing 
the sequence of the summation and the integration in Eq. (A.2) and taking into 
account the definition of Riemann's ^-function 



oo / -i \n 

E i -i L = (2 1 - s -i)C( s ), 



n-l U 



one obtains 



Cr — ioo 

It is important to note that the contours C' 2r , C' Q are the same as for K\ r \(x) after 
shifting of integration variable s — > s'. The integral (A. 3) is calculated by closing 
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the contour C to the left and summing up the residua of the integrand. For 
the contour C' 2r one should include all the poles of Y and ( -functions. But for 
positive r, < Co < 1, the pole s' — 1 of ^-function is out of the contour and 
should be excluded. 

In the case of fermion contributions the following integrals have to be calcu- 
lated: 

# > = \iC + log(^)) - § U-^m - T^-^^^n + 1), (A4) 

f?> = -^(|r'£(-?Hi-2-^') r(n r^ l/2) ((aH-»--i), (A5) 

r = 1,2..., C - Euler constant. As it is seen, in the limit uj — > (T — > oo) the 
first term in Eq. (A. 4) is dominating. Other terms give the corrections of order 

~ rpi rp2 1 

The same procedure has been applied for contributions of boson fields and for 
particular values of r the results of summations are as follows: 

± V.(«0 = ^ - ^ + f - ^ - C + 3/4) + 0((A ,A6) 

^ rr ur 2 6 16 uj 

£ i^(no;) = M _ I + "(kg 1^ - C + 1/2) + 0(u, 2 ), (A7) 
~i n c<j 2 4 uj 

V in(na;) = + a; (log — - C + 1/2) + OV), (A8) 

uj uj 

°° 1 47T 

£ #oM = -(c - log — ) + OH- (A9) 

n=l ^ u 
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